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SECTION  I 
INTRODUCTION 

Some  of  the  key  elements  In  superheterodyne  receivers,  such  as  mixers 
and  limiters,  include  semiconductor  diodes  as  their  basic  devices. 

The  nonlinear  characteristics  of  these  diodes  strongly  impact  receiver 
performance.  For  example,  insertion  loss  in  a mixer  wil1  affect 
sensitivity  while  harmonic  intermodulation  products  generated  in  the 
mixer  will  affect  the  dynamic  range  of  the  receiver;  a limiter  response 
to  multiple  signals  will  determine  the  multiple  signal  handling  capability 
of  the  receiver.  For  proper  receiver  design,  it  is  essential  to  predict 
the  performance  of  such  nonlinear  elements  under  varying  conditions. 

At  the  present  stage,  the  nonlinear  performance  of  mixers  and  limiters, 
such  as  intermodulation  distortion  and  multiple  signal  handling,  are 
determined  experimentally.  The  "spur  table"  furnished  by  mixer  manu- 
facturers is  one  example.  This  table  provides  information  at  the 
various  harmonic  frequencies  for  an  "optimal"  local  oscillator  (LO) 
drive  level.  However,  the  use  of  this  spur  table  is  quite  limited 
because  variations  in  the  LO  drive  level  will  affect  the  output  level 
of  the  intermodulation  (IMD)  products.  An  ideal  solution  to  this 
problem  would  be  to  theoretically  calculate  the  output  levels  of  the 
IMD  products  as  a function  of  the  LO  drive  level. 

In  this  report,  a mathematical  model  is  presented  by  means  of  which 
the  nonlinear  responses  of  various  diode  circuits  can  be  calculated. 

A single  diode  circuit  with  multiple  input  sources  was  simulated  to 
demonstrate  the  capabilities  of  the  model  (Section  II).  The  solution 
to  the  circuit  is  obtained  first  in  the  time  domain  using  iterative 
methods.  The  time  domain  solution  Is  then  Fourier  analyzed  to  extract 
the  frequency  components  which  are  of  Interest  to  us.  For  this  single 
diode  circuit,  single  and  two-tone  Intermodulation  distortion  levels 
were  calculated  using  this  model.  A single  ended  mixer  was  constructed 
to  verify  the  calculated  results. 

In  the  course  of  this  Investigation,  It  was  found  that  some  harmonic 
products  which  adversely  affect  superheterodyne  reception  might  in 
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some  applications  become  an  asset.  Such  Is  the  case,  for  example.  In  a 
frequency  doubling  receiver  for  detection  and  Identification  of 
pseudo-noise  biphase-shift-key  (PN-BPSK)  modulated  signals  (Appendix  A). 
Another  application  where  harmonic  Intermodulation  products  may  be 
proven  useful.  Is  the  detection  of  the  presence  of  a nonlinear  device 
within  a given  volume.  The  volume  Is  Irradiated  with  a known  signal 
and  a specified  harmonic  product  generated  by  the  nonlinear  device  is 
detected,  If  such  a device  is  present  within  that  volume. 

The  single  diode  circuit  model  was  extended  to  account  for  multiple 
diodes.  Using  the  multiple  diode  circuit  general  model,  presented  In 
Section  III-l,  we  can  simulate  a variety  of  diode  circuits  such  as 
limiters  (Section  1 1 1-2) , balanced  and  double  balanced  mixers 
(Sections  I I 1-3  and  I I 1-4  respectively).  Results  of  these  simulations 
along  with  experimental  data  are  presented  in  Section  IV. 


While  various  models  for  mixer  analysis  have  been  previously 
published  in  the  literature  (References  1 and  2),  most  of  these  models 


deal  with  the  linear  properties  of  mixers,  such  as  conversion  loss  and 


noise  figure,  and  base  their  analysis  on  a quasi-1 inear  approach  which 
assumes  a small  radio  frequency  (RF)  signal  and  a sinusoidal  local 


oscillator  waveform  on  the  mixer  diode.  These  methods  of  analysis  are 
quite  suitable  for  predicting  the  linear  performance  of  mixers  under 
small  signal  conditions  but  fail  to  predict  their  nonlinear  behavior 
under  nonidealized  conditions,  such  as  a finite  source  impedance,  large 
RF  signals,  and  diode  parasitics.  In  order  to  overcome  these  defi- 
ciencies, a more  general  model  of  a single  ended  mixer  is  presented  in 
this  section.  This  model  can  be  applied  to  predict  the  nonlinear 
behavior  of  a mixer  under  large  signal  conditions,  taking  into  account 
the  finite  source  impedance  and  the  diode  losses. 

The  single  ended  mixer  circuit  to  be  modeled  is  shown  in  Figure  la. 
The  diode  is  fed  by  a voltage  source  V.-n(t)  with  50  ft  characteristic 
impedance,  through  a circulator.  The  reflected  voltage  from  the  diode 
appears  across  the  load  impedance  1Q  at  the  output  port  of  the  circulator. 
This  particular  configuration  was  chosen  in  order  to  decouple  the 
output  from  the  input  and  ease  the  measurement  procedure  which  will  be 
later  discussed.  It  has  to  be  noted  that  this  model  does  not  include 
reactive  elements  which  ought  to  be  taken  into  account  if  proper 


representation  of  the  diode  and  package  parasitics  is  desired. 
Neglecting  reactive  elements  in  the  model  will  obviously  introduce 
discrepancies  between  the  measured  and  calculated  results.  However, 
this  model  simplifies  the  numerical  calculations  by  avoiding  the 
solution  of  nonlinear  differential  equations  which  may  include  very 
large  time  constants  to  account  for  the  d.c.  blocking  capacitors  and 
RF  chokes. 


3 


AFAL-TR-77-266 


L 


I 


The  diode  is  modeled  as  a nonlinear  conductance  in  series  with  a 


constant  resistance  (Rs).  The  current-vol tage  ( I -V)  characteristic  of 


the  diode  is  given  by: 


i = Is  [e 


a(vn-iRJ 

D S -1]  = Is(eav-1) 


(1) 


where  i is  the  current  through  the  diode,  I is  the  reverse  saturation 


current,  vD  is  the  voltage  on  its  terminals,  v is  the  voltage  on  the 


diode  junction  and  a = where  q is  the  electronic  charge,  K is  the 


Boltzmann  constant  and  n is  the  quality  factor  of  the  junction  (1  <.  n <.  2) 
and  T is  the  temperature  in  °K.  The  voltage  on  the  diode  Vp(t)  and  the 
current  through  it  i (t)  are  calculated  from  the  boundary  conditions  on 
the  diode  plane: 


vD(t)  'V'l 


i(t)  = (Vj  - vR)/Zc 


(2) 

(3) 


Vj  and  vR  are  the  incident  and  reflected  voltages  respectively  (Figure  lb). 


From  the  matched  conditions  at  the  input  port  of  the  circulator. 


VI  ■ Vtn(t>/2 


(4] 


The  relationship  between  the  voltage  on  the  diode  v^  and  the  voltage 


on  the  junction  v,  is  given  by: 


v0  = v + iR$ 


(5) 


and  the  current-vol tage  relationship  of  the  junction  itself  is  given 
by  Equation  1. 

From  Equations  1 through  5 the  following  transcendental  equation  for 
the  voltage  on  the  diode  junction  is  derived: 


F<*>  " * ♦ <Zo  + Rs>  -vin  ■ 


(6) 
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The  diode  parameters  I$,  R$  and  a are  obtained  for  dc  measurements  on 
the  diode,  and  ZQ  is  given. 

The  input  voltage  V-n(t)  is  generally  given  as  a sum  of  two  sinusoidal 
waveforms  at  the  LO  and  RF  frequencies*: 


Vin^  = VL0s1nuiL0  1 + VRFs1na)RF  1 


(7) 


Equation  (6)  is  now  solved  numerically  at  each  point  in  time  applying 
the  Newton-Raphson  method.  Let  vn  be  the  approximate  solution  of 
Equation  6 after  n iterations,  so  that: 


ctv_ 

F(,»>  ■ vn  * (Zo  * Rs>  's  <e  •»  - V1„  * E * 0 


The  solution  at  the  next  iteration  will  be  given  by: 


n+1 


= v„  - F( v )/F 1 ( v ) 


where  F'(vn)  = 


v=v. 


= l + Is“(V  Rs>  e“Vn 


(8) 

(9) 

(10) 


The  computation  stops  when  the  difference  between  the  solutions  at 
two  consecutive  iterations  falls  within  a predetermined  error  range. 

After  the  junction  voltage  v(t)  is  obtained,  the  current  through 
the  diode  i (t)  is  calculated  and  subsequently  the  voltages  Vp(t)  and 
vR(t).  The  load  voltage  is  given  by: 

VL  ’VR  = ,D  • VI  * VD  - V2  <"> 


*In  case  it  is  desired  to  calculate  the  two-tone  third-order  intermodu- 
lation (IMD)  product,  another  term  is  added  to  V^n(t): 

Vin(t)  = VL0  s1nal0  1 + VRF  (sinu)RF  * + s1nu)RF 

(7a) 

where  u^p  and  to^p  are  the  frequencies  of  the  two  tones. 
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Once  the  load  voltage  Is  known  for  an  entire  cycle,  Fourier  analysis 
Is  applied  to  obtain  the  output  voltage  at  the  various  frequencies  which 
are  of  Interest  to  us. 

2.  FREQUENCY  DOMAIN  CALCULATIONS 

The  frequencies  of  the  harmonic  intermodulation  products  generated 
In  a mixer  are  given  by: 


wmn  = Ho  + "^f! 


(12) 


where  m and  n are  positive  or  negative  integers.  All  (m,n)  combinations 
other  than  (1,  -1)  which  is  the  intermediate  frequency  (IF),  are 
considered  to  be  detrimental  to  heterodyne  reception,  but  might  be 
beneficial  in  other  applications.  The  output  amplitude  at  a particular 
IMD  product  (wmn)  Is  obtained  by  Fourier  analysis  of  vL(t): 


V-W  ■ f o'T°  t<*>  • 


-J  ui 


mn 


dt 


(13) 


where  = ^ the  period  of  vL  (t)  and  u>0  is  the  largest  common 
divisor  of  u^q  and  u^p. 


For  numerical  purposes  it  is  convenient  to  choose  the  set  (^Lo*  “rf^ 

so  as  to  minimize  TQ  in  order  to  reduce  computation  time  but  at  the 

same  time  care  has  to  be  taken  in  order  to  avoid  overlapping  of  low 

order  IMD  products,  meaning  that  u>  „ f w„  _ for  m,n  <.  6.  The 
r n1  .m^  n^.n^  — 

particular  choice  of  the  set  (u^g,wRp)  for  calculation  purposes  is 

irrelevant  to  the  actual  performance  of  the  mixer,  since  the  circuit 

model  described  in  the  previous  section  does  not  include  frequency 

dependent  elements. 

For  calculations  of  two-tone  intermodulation  levels  the  choice  of  the 
frequency  set  (u^q,  “rf*  “rF^  9u^e<*  similar  reasoning.  The  inte- 
gration (Equation  13)  is  carried  out  numerically  using  the  trapezoidal  rule. 
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The  computation  of  the  output  level  at  a given  frequency  u)mn  is  now 
complete.  From  the  given  input  voltage  ( p( t) ) the  voltage  on  the 
diode  (v(t/)  is  calculated  and  subsequently  the  voltage  on  the  load 
(vL(t))  Is  obtained  from  circuit  conditions.  The  desired  frequency 
components  of  v^(t)  are  later  calculated  by  a simple  Fourier  series 
expansion.  Results  of  calculations  using  this  model  are  presented  in 
Section  IV.  These  calculations  are  followed  by  experimental  data 
taken  on  a similar  circuit  configuration. 

In  Section  III,  the  single  diode  circuit  model  is  expanded  to 
multiple  diode  configurations. 
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SECTION  III 

MULTIPLE  DIODE  CIRCUIT  SIMULATION 

I.  GENERAL  MODEL 

Consider  the  general  circuit  configuration  shown  in  Figure  2.  The 
circuit  consists  of  a network  of  N interconnected  diodes  with  M voltage 
inputs  and  a load  impedance  (Z^).  The  diodes  may  or  may  not  be  identical 
and  the  parameters  for  each  diode  I$k,  Rgk  and  a$k  (k=l,2, N)  are 

given.  The  voltage  sources  V.  ( j=l ,2, . . . .M)  are  given  by: 

j 

V1n,  = cj  VL0  s1n(“L0  “V*  d0  VRF  Si"HF  1 * ej>  (,4> 

J 

where  c.  and  d.  are  constants. 

J J 

We  now  expand  the  single  diode  circuit  model  developed  in  Section  II 
to  solve  for  the  voltage  on  each  diode  vk(t)  and  eventually  for  the 
voltage  on  the  load  vL(t). 

The  homogeneous  loop  equation  will  have  the  form: 

F(7)  = A 7 + Z T + B Vin  = 0 (15) 

where  V is  the  vector  of  the  junction  voltages,  T is  the  vector  of  the 
diode  currents,  V.n  is  the  input  voltage  vector,  A and  Z are  N x N 
matrices,  and  B is  an  N x M matrix. 

Similar  to  Equation  9 the  approximate  solution  after  n+1  iterations 
will  be  given  by: 

(v)n+1  = (v)n  - CJ(V)ri]"7  [F(7)n]  (16) 

where  (v)n  is  the  solution  at  the  previous  Iteration  and  J(v)p  Is  the 
Jacobian  of  the  system  of  Equation  15  evaluated  at  v = (v)n 
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Figure  2.  Multiple  Diode  Circuit  Model 

The  elements  of  the  Jacobian  are  given  by. 

3Fi(v) 

j1k  = = aik  + aik  ’k 

dik  (“k  vk) 

where  i • = ^ = 0|(  I#k  e< 


(17) 

(18) 


In  summary,  once  the  matrices  A,  Z,  B,  and  the  input  voltage  vector 

V.  are  known  from  circuit  conditions,  the  solution  for  the  voltages  on 
in  _ 

the  diodes  v Is  straight-forward  using  Equation  16.  As  an  example,  the 
matrices  for  the  single  ended  mixer  case  discussed  In  the  previous 
section  are  given  by: 


A * 1 ; B * 1 ; Z * (Rs  + Z0); 
J * 1 + (R$  + Z0)  a Is  e“  v 


(19) 
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In  the  following  sections  the  above  matrices  will  be  calculated  for 
various  circuits,  namely,  balanced  and  double  balanced  mixers  and 
back-to-back  diode  limiters. 

2.  BACK-TO-BACK  DIODE  LIMITER 

The  multiple  signal  handling  capability  of  a receiver  is  strongly 
affected  by  the  performance  of  the  limiter.  The  ideal  input-output 
characteristic  of  a limiter  is  plotted  in  Figure.  3a.  Using  the  model 
developed  in  this  section  we  can  calculate  the  deviations  from  this 
ideal  characteristic  when  two  or  more  input  signals  are  present  at  the 
limiter  input. 


The  equivalent  circuit  of  a back-to-back  diode  limiter  is  shown  in 
Figure  3b.  From  loop  analysis  of  this  circuit  we  obtain  the  following 
system  of  equations: 


F1(v1,v2)  = v1  + (R$1  + Zeq)  i1  - Zeq  i2  - ^ Vin 


F2(vrv2)  = v,  + v2  + R$1  i1  + Rs2  i2  = 0 


(20) 

(21) 


where  Z_ 


Zo"V 


The  circuit  matrices  for  this  case  are  given  by: 
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3.  BALANCED  MIXER 

A simplified  circuit  of  a balanced  mixer  using  a 90°  hybrid  is  shown 
in  Figure  4a.  From  loop  analysis  of  the  equivalent  circuit  (Figure  4b) 
we  obtain  the  following  circuit  matrices: 


n ol 

-1  -( 

B -Jir/2- 

I; 

B =VF  ] 

Lo  iJ 

~T 

e-jTr/2 

1 

'WR,1 

-Z, 

Z = 


-Z, 


WRs2 


(24) 


'in 


VRF  cosujrf 


Vlo  COSa)L0  t 


The  voltage  on  the  load  is  given  by: 


vL  = (1,  - 12)  ZL 


(25) 


Note  that  the  input  voltage  vector  V^n  is  chosen  in  such  a manner 
as  to  enable  rapid  convergence  of  the  solution  at  t=0.  Rapid  convergence 
at  t=0  Is  assured  when  the  diodes  are  back  biased  or  at  a very  low 
forward  conduction  state. 

For  a balanced  mixer  using  a 180°  hybrid,  the  A and  Z matrices  are 
the  same  as  in  the  90°  hybrid  case."  The  B matrix  and  the  Input  voltage 
vector  V^n  will  be  given  in  this  case  by: 


B - VF_ 
B - T 


-1  -1 

m m 

VL0  s1nu*L0  t 

1 

* Vin  = 

-1  -1 

v 

RF  slncopp  t 

(26) 
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Vin2«*> 


Vjn = VRF  cos  WRF  t 
vin2<*)  = VL0  cos  WL0  f 


(o)  Simplified  Circuit 


Figure  4.  90°  Hybrid  Balanced  Mixer 


AFAL-TR-77-266 


I 


4.  DOUBLE  BALANCED  MIXER 

The  simplified  circuit  of  a double  balanced  mixer  is  shown  in 
Figure  5a.  This  mixer  is  analyzed  using  the  equivalent  circuit  in 
Figure  5b.  A and  Z are  in  this  case,  4x4  matrices  and  B is  a 
4x2  matrix. 


"l 

1 

0 

o" 

”-l 

o" 

0 

1 

1 

0 

0 

-1 

0 

1 

0 

0 

; B = 

-*5 

-*s 

_1 

1 

1 

1 _ 

. 0 

0_ 

Figure  5.  Double  Balanced  Mixer 
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SECTION  IV 

COMPUTER  EXPERIMENTS  ON  DIODE  CIRCUITS 

1.  CALCULATIONS  AND  MEASUREMENTS  ON  A SINGLE  ENDED  MIXER 

Two  specific  applications  for  the  single  ended  mixer  model  were 
investigated  in  particular.  The  first  case  is  the  detection  of  the 
presence  of  a nonlinear  device  within  a given  volume  by  irradiating 
this  volume  with  a known  signal  and  detecting  a specified  harmonic 
product  of  the  transmitted  signal. 

This  method  is  widely  applied  in  commercial  antishoplifting  systems 
where  a diode  is  attached  to  the  protected  merchandise  and  irradiated 
by  a signal  at  f . A receiver  at  the  entrance  to  the  store  detects  the 
second  harmonic  at  2f  generated  by  the  diode  and  activates  an  alarm. 

A variation  of  this  method  can  be  applied  when  it  is  preferred  to  have 
the  harmonic  product  generated  by  the  nonlinear  device  in  the  same 
frequency  range  as  the  transmitted  signal.  For  this  purpose,  two  equal 
amplitude  signals  at  f^  and  f2  are  transmitted  so  that 

V i n ( t ) = VQ  (sin  u^t  + sin  u^t)  (28) 

The  product  | 2f 1 - f2|  (or  | 2f 2 - f 1 | ) is  then  detected  if  a nonlinear 
device  is  present  within  the  irradiated  volume. 

The  model  developed  in  Section  II  was  applied  to  evaluate  the  efficacy 
of  this  method.  The  calculated  and  measured  power  output  at  various 
intermodulation  products  as  a function  of  the  power  input  ( P ^ n ) is 
shown  in  Figures  6 and  7 respectively. 

The  experimental  results  were  taken  on  a homemade  mixer  constructed 
by  mounting  a detector  diode  in  a shorted  N-type  connector  without  the 
benefit  of  any  matching  networks.  The  experimental  setup  used  is  shown 
In  Figure  8.  The  static  I-V  characteristics  of  the  diode  were  measured 
on  a curve  tracer  and  the  diode  parameters  a,  Is»  and  R$  were  matched 
to  this  curve  and  used  in  the  computer  simulation. 
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Pin  (dBm) 


Figure  6 


Calculated  Power  Output  (Pout)  of  (m,n)  IMD  Products  as 
a Function  of  the  Power  Input  (Pin).  (pi=P2=pin^ 


igure  8.  Experimental  Setup  for  the  Measurement  of  the  Power  Output 
of  IMD  Products  Generated  by  the  Mixer  (P.=P?) 
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The  results  of  this  simulation  (Figure  6)  show  that  for  some  levels 
of  input  power,  some  of  the  IMD  products  completely  disappear.  A similar 
phenomenon  was  observed  in  the  experimental  data,  as  shown  In  Figure  7. 
Although  in  this  case  the  output  levels  do  not  drop  as  drastically  as  in 
the  simulated  case,  dips  in  the  power  output  of  the  various  IMD  products 
were  seen  to  occur  for  given  input  power  levels.  The  Numerical  dis- 
crepancies between  the  calculated  and  measured  results  could  be  attributed 
to  the  fact  that  the  computer  model  does  not  include  frequency  effects, 
while  the  experimental  setup  is  frequency  limited.  In  addition,  the  I-V 
diode  characteristic  is  not  strictly  exponential  as  assumed  in  the 
computer  model. 

Another  possible  application  for  this  model  is  to  predict  the  level 
of  the  second  harmonic  generated  in  a diode  in  order  to  evaluate  the 
feasibility  of  detection  and  identification  of  pseudo-noise  BPSK 
modulated  signals  by  a frequency  doubling  technique.  The  details  of 
this  investigation  are  presented  in  Appendix  A. 

The  two-tone  third-order  intercept  point  in  a mixer  can  be  predicted 
as  a function  of  the  local  oscillator  drive  level  using  the  model 

presented  here.  For  example,  the  intercept  point  for  the  single- 

ended  mixer  previously  simulated  was  calculated  for  three  values  of 
LO  drive  level.  The  results,  along  with  the  RF  to  IF  conversion  loss  at 

these  LO  levels,  are  shown  in  Table  1. 

TABLE  1 


TWO-TONE  THIRD  ORDER  INTERCEPT  POINT  AND  INSERTION  LOSS  AT 
VARIOUS  LO  DRIVE  LEVELS.  (SINGLE-ENDED  MIXER) 


^dBm) 

THIRD  ORDER  INTERCEPT  (dBm) 

INSERTION  LOSS  (dB) 

3 

11.1 

11.3 

6 

14.5 

10.9 

9 

17.7 

10.7 
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It  Is  seen  that  In  this  range,  an  Increase  of  3 dB  in  the  LO  drive 
level  will  increase  the  intercept  point  by  more  than  3 dB  while  slightly 
improving  the  Insertion  loss  of  the  mixer.  Thus,  by  using  this  model, 
an  LO  power  can  be  calculated  which  would  be  optimal  for  our  purposes, 
be  it  maximum  linearity,  minimum  conversion  loss,  suppression  or 
enhancement  of  a particular  IMD  product,  etc. 

2.  BACK-TO-BACK  DIODE  LIMITER 


The  performance  of  a back-to-back  diode  limiter  with  an  input  con- 
sisting of  two  signals  at  f^  and  f^,  was  simulated  using  the  formulation 
presented  in  Section  1 1 1-2.  The  output  power  at  f^  (PQut  ) was  calcuated 

as  a function  of  the  input  power  at  f,  (P.  ) with  the  input  power  at 

i 1 ri-j 

f2  ^Pin  ^ as  a Parameter-  Results  of  this  calculation  are  shown  in 

Figure  9.  The  simulated  diodes  were  considered  lossless  (R$=0),  but  for 
one  case  where  a spurious  series  resistance  was  added  (R$  = 5ft)  in  order 
to  evaluate  its  effect  on  the  performance  of  the  limiter  (dotted  line). 


As  seen  in  Figure  9,  the  presence  of  a series  resistance  in  the  diode 
is  detrimental  to  the  limiting  effect.  The  increase  in  current  through 
the  diodes  due  to  increasing  input  power  yields  a higher  voltage  on  the 
diode  terminals,  canceling  the  limiting  effect  of  the  junction  itself. 


Figure  9 also  shows  the  effect  of  the  presence  of  a second  signal  at 


the  input  of  the  limiter.  The  attenuation  is  linear  with  P 


in,  P,n2  * 


P„„.  where  Pm,w  is  the  maximum  output  power 
max  max 


This  is  the  "capture  effect" 
of  the  limiter,  where  the  stronger  signal  prevails  and  the  weaker  signal 
is  further  attenuated.  This  effect  is  clearly  seen  In  Figure  10,  where 
the  output  power  of  both  signals  is  plotted. 
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3.  90°  HYBRID  BALANCED  MIXER 

A simulation  of  a 90°  hybrid  balanced  mixer  with  identical  diodes 
was  carried  on  using  the  equations  in  Section  II 1-3.  Two  cases  were 
investigated: 

a.  Insertion  loss  and  two-tone  third-order  intercept  as  a function 
of  the  local  oscillator  power  level  (Figure  11). 

b.  IF  and  single-tone  IMD  output  levels  as  a function  of  RF  input 
power  (Figure  12). 

As  in  the  single-ended  mixer  case,  increasing  the  LO  power  level 
will  greatly  improve  linearity  (higher  intercept  point)  while  slightly 
reducing  insertion  loss  (Figure  11).  The  IF  and  IMD  plots  (Figure  12) 
are  similar  to  those  supplied  by  mixer  manufacturers  but  for  the  "dips" 
in  the  output  levels  of  the  2x2  and  1 x 2 products,  (m  x n = |mfLQ  - nfRp|). 
These  dips  appear  in  this  simulation  at  high  input  power  levels  generally 
not  measured  by  the  manufacturers.  No  physical  explanation  for  this 
phenomenon  has  been  found  by  the  authors.  Mathematically,  the  load 
voltages  at  these  products,  which  are  calculated  from  a Fourier  series 
expansion,  change  sign  at  these  input  RF  levels. 

Additional  examples  of  applications  of  the  model  are  presented  in 
Appendix  B,  along  with  program  usage  Instructions  and  the  program  listing. 
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INSERTION  LOSS  (dB) 


LO  POWER  (dBm) 


Figure  11.  Calculated  Insertion  Loss  and  Two-Tone  Third-Order  IMD 
Intercept  Point  as  a Function  of  the  LO  Power,  for  a 90° 


TWO-TONE  3rd  ORDER  INTERCEPT  POINT  (dBm) 


SECTION  V 


CONCLUSIONS 

The  multi-diode  circuit  model  for  computer  simulation  developed 
during  the  course  of  this  Investigation  has  proved  to  be  useful  for 
qualitative  analysis  of  mixers,  as  well  as  other  diode  circuits. 
Investigation  of  possible  applications  of  mixers  other  than  heterodyne 
mixing  can  be  easily  carried  out  using  this  model. 

An  effort  should  be  made  to  improve  the  model  in  order  to  achieve 
quantitative  agreement  between  measured  and  calculated  results.  This 
effort  should  include  diode  parasitics  to  obtain  more  realistic  results. 
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PSEUDO-NOISE  BPSK  SIGNAL  DETECTION  BY  HETERODYNE  MIXING 
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Detection  of  the  presence  of  Pseudo-Noise 
Bi-Phase  Shift  Keyed  signals  rom  uncooperative 
sources  is  extremely  difficult  due  to  their 
noise  like  and  suppressed  carrier  characteris- 
tics. 

A method  which  facilitates  detection  of 
such  signals  Is  presented.  The  method  consists 
of  squaring  the  incoming  signal  by  means  of  a 
non-linear  device  in  the  presence  of  a pumping 
signal.  The  squaring  process  eliminates  the  bi- 
phase modulation,  thus  obtaining  a CW  signal  at 
the  second  harmonic  of  the  carrier  frequency. 

By  pumolpg  the  non-linear  device  with  a local 
oscillator  signal,  significant  enhancement  of 
the  second  harmonic  output  is  obtained,  with  a 
distinctive  peak  in  output  power  at  a given  jxmp 
level.  Moreover,  this  method  makes  it  possible 
to  down-convert  tne  second  harmonic  signal  to  a 
lower  frequency  without  considerable  loss, 
faci litating  the  post-processing. 

The  described  method  was  tested  success- 
fully both  by  computer  simulation  and  laboratory 

experiments. 


INTRODUCTION 

The  non-linear  properties  of  a heterodyne 
mixer,  which  adversely  affect  reception  by 
generation  of  spurious  signals,  might  in  some 
applications  become  an  asset.  Such  Is  the  case 
presented  here  in  which  the  mixer  non-linearities 
ars  Investigated  as  a means  to  the  detection  of 
the  presence  of  Pseudo-Noise-BPSK  signals  from 
unr.ooptratlve  sources. 

The  time  function  of  a pseudo-nolse-bi- 
phase  modulated  signal  is  given  by: 


This  spectrum  is  obviously  similar  to  that  of 
a pulsed  RF  signal  with  a small  duty  cycle.  An 
operator  monitoring  the  spectrum  and  delecting  a 
signal  with  such  a power  distribution  on  his  dis- 
play could  not  immediately  decide  if  it  is  a PN  or 
a pulse  modulated  signal.  If  time  is  no  factor, 
however,  a FN  signal  could  be  identified  by  its 
suppressed  carrier  characteristics.  Fig  1(b). 

An  additional  method  for  identification  of  a 
PN  signal  as  such  is  suggested  here.  The  method 
consists  of  obtaining  a CW  signal  at  the  second 
harmonic  of  the  carrier  frequency  by  squaring  the 
incoming  signal  in  the  presence  of  a pump  (local 
oscillator).  Applying  this  method,  it  is  also 
possible  to  obtain  a down-converted  CW  signal  by 
heterodyne  mixing  of  the  pump  frequency  with  the 
second  harmonic  of  the  carrier. 

MATHEMATICAL  FORMULATION 

The  squaring  of  the  incoming  signal  is  accom- 
plished by  means  of  a n.Qp-llnear  device,  such  as  a 
semi-conductor  diode.  The  output  from  such  a 
device  can  be  represented,  for  a small  signal 
regime,  as  e power  series  of  the  form: 

Vout(t>  ■ * \ Vi>>  <2> 

k*o 

Expanding  the  quadratic  terra  with  V (t)  as  given 
in  (1),  we  get:  n 

Vin2(t)  - Vc2  cos2  {V  + ♦<*»  ‘ 

» ~ Vc2  (cos  £2^1  + 2A(t)]  4-1) 

..  .but  since  p(t)  is  either  ci  or  *,  therefore 


V(t)  • V£  cos  (a£t  4 »(t))  (1) 

vhera  u Is  the  carrier  frequency  and  the  phase 
♦ (t)  chSngcs  In  a pseudo  random  fashion  between 
0 and  *,  controlled  by  a direct  sequency  modu- 
lator. 

The  powar  spcctrun  of  such  a signal  which 
appears  to  be  noise  with  a (sinx/x)^  distribution 
is  shown  In  Fig  1(a). 


Vtn 


r Vc2  [co.2*ct4l] 


(3) 


Thus  eliminating  the  phase  modulation  of  the  orig- 
inal signal  and  obtaining  a CW  signal  at  2w  with 
an  amplitude  of  A?V  /2;  if  we  take  into  account 
the  fourth  order  term  also,  the  amplitude  at  the 
second  harmonic  will  be 
2 
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Fig  Is  Pow«r  Spectrum.  of  a PN-5PSK  signal 

(10  dB/dlv) . f * 2200MHI,  1 MHz  clock,  127  bits 

sequence  length. 


Fig  2:  Calculated  power  output  at  2o>  -u  as  a 
function  of  the  pump  power  P , with  tfie  carrier 
power  Pc  as  a parameter.  " 
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Fig  *1  Measured  power  output  at  2u  -w  at  a 
function  of  the  pump  power  P , with^the  carrier 
power  P£  as  a parameter.  p 
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This  process  is  bound  to  be  lossy  due  to  its  quad- 
ratic characteristic,  especially  at  snail  signal 
levels,  where  this  model  applies.  Some  improve- 
ment can  be  obtained,  however,  by  sviperiniposing 
or.  the  Incoming  signal  a pump  vita  an  amplitude 
of  V at  a frequency  of  c . so  that  the  input 
voltEge  to  the  nen-linear*'device  will  be  now: 

V.  (t)  • V cosjw  t ♦ o(t)*l  4 V cos  w t (4) 
n c w c p p 

Let  us  now  calculate  the  second  and  fourth  order 
terns  in  (2)  and  extract  tne  second  harnonic 
component : 

A,V.  2 » A^{V  cos  fu)  t + <(t)l  + V cosj  t}fc  - 
win  «.  c c J PP 

A,V  2 

• ~,C  cos2u  t + •••  (5) 

i c 


L a 4 \k  2 2 

A. V -C4  V + J 4 V ) cos2u-  t +• 
4 in  c -y-  c p c 


Adding  (5)  and  (6)  we  get: 

V <2u  ) - y-  V 2(1  + -*•  V 2 + — A V 2)  (7) 

out  c 2 c 'A^  c A9  p 


Thus,  an  additional  term  is  added  to  the  amplitude 
of  the  pump  signal  V . If  A,  and  have  the 
same  sign,  the  anpliEude  of  the  CV  signal  at  2u^ 
will  increase  with  the  pump  level.  Fortunately, 
as  will  be  seen  later,  this  is  indeed  the  case  up 
to  some  level  of  V where  the  small  signal 
approximation  collapses. 


From  the  power  scries  in  (2),  we  can  also 
obtain  a down  converted  CW  signal  at  a frequency 
w.  - ■ ju.  - 2u»  i , by  calcinating  tne  third  order 
m*  P ^ 

term  in  the  series.  The  amplitude  of  this  signal 

will  be: 

V „ (U.  ,)  - ^A.V  2V  (8) 

out  1,2  23c  p 


It  has  to  be  noted  though,  that  a CW  signal 
of  amplitude  V at  the  carrier  frequency  w will 
generate  a second  harmonic  of  the  sar.e  amplitude 
as  a PN*LPSK  signal,  therefore,  in  order  to  posi- 
tively identify  such  a signal,  the  carrier  band  has 
to  be  monltcrcd  simultaneously. 


The  aoove  theoretical  results,  which  inJicate 
that  an  improvement  in  the  second  harnonic  output 
can  be  achieved  in  the  presence  of  a pumping 
signal  were  confirmed  both  by  numerical  calcula- 
tions and  by  laboratory  measurements. 


NUMERICAL  RESULTS 


In  order  to  confirm  the  above  theoretical 
results,  a computer  simulation  of  a single-ended 
mixer  was  carried  out.  The  model  of  the  mixer 
consists  of  a semiconductor  diode  with  an 
exponential  I-V  characteristic  and  a series 
resistance  R#  representing  the  diode  losses.  The 


diode  is  fed  by  a 50“  signal  source  through  a 
circulator  in  order  to  isclate  the  input  from  the 
output.  The  DC  parameters  of  a detector  diode, 
which  was  later  used  in  our  experiments,  were 
obtained  from  the  measured  1-V  curve.  From  tills 
model  a transcendental  equation  for  the  voltage  on 
the  load  as  a function  of  tine,  is  obtained.  This 
equation  is  solved  using  the  Nevton-Raphson 
method.  After  the  time  function  is  calculated, 
Fourier  analysis  is  applied  in  order  to  calculate 
the  various  spectral  components  of  interest. 

In  the  case  reported  here,  the  power  output 
at  the  second  harmonic  2utc  and  at  the  down  con- 
verted frequency  u.  - ■ j 2u.  -u.  was  calculated  as 

1 , A.  C p 

a function  of  P , the  pump  power  level,  with  P , 
the  carrier  power  level,  as  a parameter.  The 
results  are  shown  in  Fig  2.  It  Is  observed  that 
the  outputs  at  the  second  harmonic  and  at  the  down 
converted  product  exhibit  a peak  at  a certain  pump 
power  level,  decreasing  for  higher  pump  levels, 
leading  to  the  conclusion  that  the  small  signal 
power  series  model  developed  in  the  previous 
section  is  not  valid  anymore  at  these  input  levels. 
At  the  peak,  the  second  harnonic  output  is  about 
34  dB  over  the  uon-pumped  level  in  this  particular 
simulation,  a significant  improvement. 

EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 

Using  the  experimental  setup  shown  in  Fig  3, 
measurements  of  the  power  output  at  the  second 
harmonic  and  at  the  down  converted  product  w ^ 
were  carried  on,  as  a function  of  the  pump  * 
power  level.  These  results  are  shown  in  Fig  4. 
These  plots  show  a qualitative  similarity  to  these 
obtained  by  computer  simulation,  alas,  numerically 
they  differ.  Especially*  the  boost  in  second 
harmonic  generation  under  optimal  pump  power  is 
now  only  18  dB  above  the  non-pumped  level.  These 
differences  might  be  attributed  tc  the  crudity  of 
our  computer  model  which  does  not  take  into 
account  frequency  effects;  thus,  such  discrepancies 
might  well  be  expected  at  microwave  frequencies. 
However,  the  experimental  results  corroborated  the 
computer  calculations,  at  least  qualitatively, 
showing  the  feasibility  of  this  method  for  positive 
identification  of  FN-BPSK  signals,  provided  that 
simultaneous  monitoring  is  carried  on  at  the 
carrier  band. 
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APPENDIX  B 

COMPUTER  PROGRAM  USAGE 


The  computer  program  MIXER  is  written  in  FORTRAN  IV  and  intended  for 


interactive  use. 

Data  is  entered  from  the  terminal  as  answers  to  queries  by  the 
program. 

1 . INPUT 

MIXER  accomodates  the  following  options  according  to  the  specified 
number  of  diodes  in  the  circuit,  ND: 

ND=1  : Single  ended  mixer. 

ND=2  : Limiter  or  balanced  mixer. 

ND=4  : Double  balanced  mixer. 

If  ND=2  the  program  asks  whether  the  circuit  to  be  analyzed  is  a 
balanced  mixer  or  a limiter.  If  it  is  a mixer,  the  type  of  hybrid  used 
(90  or  180)  has  to  be' specified. 

For  ND=2  the  option  exists  to  specify  non-identical  diodes  in  the 
circuit.  This  feature  is  useful  to  check  for  effects  of  unbalance  on 
supposedly  balanced  configurations. 

In  case  of  trouble,  it  is  desirable  to  obtain  the  output  in  the  time 
domain  in  order  to  check  for  nonconvergence  of  the  solutions.  This 
option  is  available  by  answering  YES  to  the  "TIME  DOMAIN  ONLY?"  query. 

In  this  case,  the  voltages  on  the  diodes  and  the  convergence  error  at 
each  Iteration  could  be  obtained  by  answering  YES  to  the  "INTERMEDIATE 
ITERATION  OUTPUT"  query. 

For  frequency  domain  simulations  the  option  Is  between  calculation  of 
the  one-tone  harmonic  intermodulatlon  products  (SPURS)  or  the  two-tone 
third-order  Intermodulatlon  product  (IMD).  The  levels  of  the 
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intermodulation  products  can  be  calculated  as  a function  of  the  LO  level 
(PLO)  with  the  signal  level  (PSIG)  as  a parameter  (1),  or  as  a function 
of  PSIG  with  PLO  as  a parameter  (2),  or  as  a function  of  PIN  (3),  where 
PIN=PLO=PSIG,  for  calculation  of  two-tone  sideband  levels. 

After  all  the  options  for  the  run  have  been  entered,  the  data  is 
specified  next  in  the  following  order: 

Diode  parameters:  ALPHA  (1/VOLTS),  saturation  current  (IS)  (AMPERES), 
series  resistance  (RS)  (OHMS). 

Input  and  output  impedances  (ZO  and  ZL)  (OHMS). 

LO  and  RF  frequencies  (FLO  and  FSIG)  (HZ). 

If  IMD  was  specified,  enter  FSIG2  (second  tone  frequency). 

Power  Levels  and  power  steps  (dBm): 

Option  (1):  PSIGMIN,  PSIGMAX,  STEPSIG,  PLOMIN,  PLOMAX,  STEPLO. 

Option  (2):  PLOMIN,  PLOMAX,  STEPLO,  PSIGMIN,  PSIGMAX,  STEPSIG. 

Option  (3):  PMIN,  PMAX,  PSTEP. 

(Example:  Calculate  SPURS  as  a function  of  PSIG  for  -30  <.  PSIG 
£.  -20  dBm  each  2 dB  while  PLO  varies  from  0 to  10  dBm  in  5 dB  steps: 
PLOMIN,  PLOMAX,  STEPLO,  PSIGMIN,  PSIGMAX,  STEPSIG  = 0,  10,  5,  -30,  -20,  2. 
For  time  domain  only  calculations  only  one  set  of  power  levels  is 
specified:  P10,  PSIG. 

If  "SPURS"  was  specified  the  number  of  Spurs  desired  (NPRODS)  is 
entered  next.  The  order  (m,n)  of  the  spurs  is  entered  as  follows: 
ml,  nl  .m2,  n2...  m is  the  coefficient  of  FLO  while  n is  the  coefficient 
of  FSIG.  For  example,  the  IF  will  be  given  by  (1,  -1);  the  LO  feed- 
through is  given  by  (1,0)  and  the  RF  feedthrough  is  given  by  (0,1). 

For  time  domain  only  simulation,  the  time  limit  of  the  calculation 
TMAX  has  to  be  specified  in  seconds.  For  frequency  domain  calculations 
TMAX  is  internally  calculated. 

Next,  the  integration  step  DELTAT  is  entered  (seconds). 

The  maximum  number  of  iterations  and  the  maximum  allowed  convergence 
error  are  entered  last. 
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2.  INPUT  CONSIDERATIONS 

a.  Choice  of  Frequencies 

As  explained  in  the  text,  the  choice  of  the  set  (FLO,  FSIG)  does  not 
affect  the  actual  performance  of  the  mixer  but  will  affect  the  results 
of  the  simulation  if  not  carefully  chosen.  The  sampling  windown  TMAX 
is  internally  calculated  as:  TMAX  = 1 / | FLO-FSIG | for  spur  calculations. 
In  order  to  avoid  erroneous  results  1/TMAX  has  to  be  a common  divisor 
of  FLO,  FSIG  and  all  their  byproducts.  A reasonable  choice  for  the  set 
(FLO,  FSIG)  would  be  (20,  19)  so  that  TMAX=1 . This  set  also  assures 
that  there  will  be  no  overlapping  of  low-order  intermodulation  products. 

For  two-tone  third-order  intermodulation  calculations,  the  choice  of 
the  set  (FLO,  FSIG,  FSIG2)  is  more  critical  yet  since  there  are  more 
intermodulation  products  to  contend  with  and  overlapping  should  be 
avoided.  TMAX  in  this  case  is  calculated  as  1/ | FSIG-FSIG2 | and  has  to 
be  a common  divisor  of  FLO,  FSIG,  FSIG2  and  all  their  byproducts.  Again 
it  is  convenient  to  choose  TMAX=1  as  before.  We  have  found  that  the 
set  (30,  12,  11)  works  reasonably  well. 

b.  Time  Increment  DELTAT 

According  to  the  Sampling  theorem,  the  sampling  rate  has  to  be  at 
least  twice  the  maximum  frequency  present  in  the  spectrum  of  the  analyzed 

_3 

signal.  It  was  found  that  for  the  above  frequency  choices,  DELTAT=10  s 

-4 

for  spur  calculations  and  DELTAT=10  s for  two-tone  third-order  IMD 
calculations  give  accurate  results.  It  is  advisable  to  check  doubtful 
results  by  running  the  program  again  with  the  same  data  and  a reduced 
DELTAT. 


c.  Number  of  Iterations  and  Maximum  Error 

Generally  the  solutions  converge  in  less  than  five  iterations  with  an 
error  of  less  than  10‘7.  It  is  advisable  to  specify  the  number  of 
iterations  as  10  or  greater  and  the  maximum  allowable  error  as  10~7  or 
less. 
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3.  OUTPUT 


The  frequency  domain  output  (spur  levels)  is  filed  under  TAPF1  while 
the  time  domain  output  is  filed  under  TAPE2. 

j 

Up  to  six  IMD  products  can  be  accomodated  on  an  iteractive  terminal. 
If  the  number  of  spurs  (NPRODS)  is  greater  than  six,  the  output  has  to 
be  diverted  to  a line  printer.  Up  to  10  spurs  can  be  requested. 

4.  CONTROL 

Mixer  uses  the  matrix  inversion  subroutine  LINV1F  which  is  part  of 
the  IMSL  library.  To  attach  this  library  the  following  command  is  used: 
ATTACH,  IMSL,  ID=X654321,  SN=AFIT. 

To  run  the  program:  XEQ,  LDSET,  LIB=IMSL,  LOAD-LFN  where  LFN  is 
the  file  name  under  which  the  object  deck  of  MIXER  is  filed.  If 
immediate  output  on  the  terminal  is  desired  use:  DISCONT,  TAPE1  (or 
TAPE2)  before  the  XEQ  command. 

5.  EXAMPLES 


In  the  following  pages  various  examples  of  the  use  of  this  program 
are  presented. 
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■»*  * . I jppn  COPY  FUK8ISHED  10  M>C 

• i • • 


PjUkCIXCABa 


EXAMPLE  1 

CALCULATION  OF  TWO-TONE  THIRD-ORDER  IMD  LEVELS 
IN  A SINGLE-ENDED  MIXER 
(a)  Input 


**t  NON-LINEAR  CIRCUIT  SIMULATION  *** 

ENTER  OPTIONS  FOR  THIS  RUN: 

ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT. 

(TYPE  1.2.  OR  4)  ND-_1_ 

TIME  DOMAIN  ONLY?  (YES  OR  NO):  N£ 

CALCULATE  SPURS  OR  2-TONE  3RD  ORDER  IMDT 
TYPE  SPURS  OR  IMD : IMP 

DETERMINE  WHAT  KIND  OF  CALCULATION  DESIRED. 

(l)POUT-F(PLO) . (2)P0UT-F(PSIG).  (3>P0UT-F(PIN) .PIN-PLO-PSIG. 
TYPE  I .2  OR  3. : X 

ENTER  RUN  DATA: 

ENTER  DIODE  «1(  ALPHA. IS. RS-  40. . 1 . E-12. 10. 

ENTER  INPUT  AND  OUTPUT  IMPEDANCES  ZO.ZL-  50. .50. 

ENTER  FREQUENCIES  FLO.FSIG  (HZ)l  30. .11. 

ENTER  2ND  TONE  FREQUENCY  FSIG2:  12. 

ENTER  POWER  LEVELS  (IN  DBM ) : 

PSIGMIN.PSIGMAX.STEPSIG.PLOMIN.PLOMAX.STEPLO:  -20.-10.10.3.9.3 

ENTER  SAMPLING  INTERVAL  DELTAT  (SECS):  l.E-4 

— 

MAX.  NUMBER  OF  ITERATIONS.  MAX.  ALLOWED  ERROR- 10 .l.E-7 

***  END  OF  DATA  INPUT. BEGIN  EXECUTION.  ««* 

***  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE1 • *** 

***  FIND  TIME  DOMAIN  OUTPUT  ON  TAPE2.  *** 

| I 


SIMULATION  OF  A SINGLE  ENDED  MIXER 
•CALCULATE  THE  POWER  OUTPUT  OF  2-TONE  3RD  ORDER  IMD  PRODUCTS 
AS  A FUNCTION  OF  PLO 

DIODE# 1 • ALPHA-  #0.00  IS-  I.OOE-12  RS-  10.00 
FREQUENCIES!  FLO-  3.00E+01  FSIG-  l.lOE+Ol  FSIQ2-  1 • 20E+01  HZ 
SOURCE  IMPEDANCE-  50.00  LOAD  IMPEDANCE-  50.00 

RUN  PARAMETERS:  DELTAT-  1.00E-04  NITER-10  ERMAX-  1.0E-07 


*********************************************************************** 


PSIG-  -20.00  DBM 


* PLO 

* IMD  * 

IF 

« 

* 3.00 

* -84.19  * 

-28.27 

* 

* 6.00 

* -92.93  * 

-26.86 

» 

» 9.00 

« -100.24  * 

-26.22 

* 

***««*«*««*#«*«««***#»******#*««****t*«*tM***«»*t*«t***t«*********«*** 

PSIG-  -10.00  DBM 

» PLO 

* IMD  « 

IF 

» 

* 3.00 

* -53.14  * 

-18.69 

* 

* 6.00 

* -62.61  « 

-16.98 

» 

* 9.00 

* -70.14  * 

-16.27 

* 

**«r*4E****««t«c4c#c******#c****«*««c**«c*****«[«M(c**4E*4i*4t«***4[**«4c«)(c*«;4r»«[4t******* 


WIS  rAQB  IS  BEST  QUALITY  PRACTICABLE 
AFAL-TR-77-266  ' ' ltati^OOPY  FURNISHED  TO  DOC  

EXAMPLE  2 

CALCULATION  OF  SIDEBAND  GENERATION 
IN  A SINGLE-ENOED  MIXER 
(a)  Input 


**«  NON-LINEAR  CIRCUIT  SIMULATION  *** 

ENTER  OPTIONS  FOR  THIS  RUN? 

ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT. 

(TYPE  It 2*  OR  4)  ND=1_ 

TIME  DOMAIN  ONLY?  (YES  OR  NO)>  NO 

CALCULATE  SPURS  OR  2-TONE  3RD  ORDER  IMDT 
TYPE  SPURS  OR  IHDt  SPURS 

DETERMINE  WHAT  KIND  OF  CALCULATION  DESIRED. 

( 1 >POUT“F (PLO) • (2)POUT«F(PSIG).  (3>P0UT=F(PIN> » PIN-PLO-PSIG. 
TYPE  1 .2  OR  3.(  3 

ENTER  RUN  DATA? 

ENTER  DIODE  *11  ALPHA. IS.RS=  40. . 1 .E-12. 10. 

ENTER  INPUT  AND  OUTPUT  IMPEDANCES  ZO.ZL-  50. .50. 

ENTER  FREQUENCIES  FLO.FSIG  <HZ> 1 20. .19 

ENTER  POWER  LEVELS-  (IN  DBM) l 
PMIN.PMAX.PSTEP— 30.0.3 

ENTER  NUMBER  OF  SPURS  DESIRED.  (1  TO  10)  J 4_ 

ENTER  SPUR  ORDER  M.N.  (M*FLO+N*FSIG> 1 1.-2  2.-3  3.-4  4.-5 

ENTER  SAMPLING  INTERVAL  DELTAT  (SECS) I l.E-3 

NAX.  NUMBER  OF  ITERATIONS.  MAX.  ALLOWED  ERROR -10.1. E-7 

***  END  OF  DATA  INPUT. BEGIN  EXECUTION.  *** 

***  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE 1 • *** 

***  FIND  TIME  DOMAIN  OUTPUT  ON  TAPE2.  *** 
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t.  jWtSLPAQB  IS  BEST  QUALITY  PRACTICABLE 
ISIOif  COPY  FURNISHED  TO  UDC 


EXAMPLE  2 


(b)  Output 


SIMULATION  OF  A SINGLE  ENDED  MIXER 
♦CALCULATE  THE  POWER  OUTPUT  OF  SPURIOUS  PRODUCTS 
AS  A FUNCTION  OF  PIN.PIN-PSIG-PLO. 


HZ 


DIODE411  ALPHA-  40.00  IS-  1.00E-I2  RS-  10.00 
FREQUENCIES:  FLO-  2.00E+01  FSIG-  1.90E+01  FSIG2-  0. 

SOURCE  IMPEDANCE-  50.00  LOAD  IMPEDANCE-  50.00 

SPURIOUS  5 < 1.-2X  2.-3H  3.-4H  4.-5) 

RUN  PARAMETERS:  DELTAT-  1.00E-03  NITER-10  ERMAX-  1.0E-07 


**iM************t****************************»**«********************** 


PIN-PLO1 
PIN  * 
-30.00 
-27.00 
-24.00 
-21.00 
-18.00 
-15.00 
-12.00 
-9.00 
-6.00 
-3.00 
0.00 


'PSIG 
< 1.-2) 
-218.75 
-208.63 
-197.47 
-184.36 
-167.85 
-145.74 
-114.97 
-72.07 
-37.61 
-25.72 
-20.48 


< 2.-3) 
-250.68 
-234.80 
-218.22 
-200.20 
-179.47 
-154.02 
-120.79 
-76.32 
-45.79 
-53.35 
-34.58 


< 3.-4) 
-281.83 
-266.44 
-244.41 
-221.04 
-195.44 
-165.83 
-129.28 
-82.69 
-69.54 
-45.08 
-54.56 


< 4.-5) 
-277.74 
-275.09 
-281.21 
-246.47 
-215.04 
-180.72 
-140.25 
-91.30 
-63.08 
-59.07 
-45.70 


**«********»**««««««*«**«««»«»»««*«***««***«*«»***«*«*******«****»***«* 


. 

I 


I 
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EXAMPLE  3 

CALCULATION  OF  ONE-TONE  IMD  LEVELS  IN  A 90°-HYBRI0 
BALANCED  MIXER  AS  A FUNCTION  OF  PLO 
(a)  Input 

*»*  NON-LINEAR  CIRCUIT  SIMULATION  *** 

ENTER  OPTIONS  FOR  THIS  RUN! 

ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT. 

(TYPE  If 2.  OR  4)  ND»2 

MIXER  OR  LIMITERT  (TYPE  N OR  Lit  M 

90  OR  180  DEGREE  HYBRID?  TYPE  VO  OR  1801  90 

IDENTICAL  DIODES?  (YES  OR  N0>!  YES 

TIME  DOMAIN  ONLY?  (YES  OR  NO) I NO 

CALCULATE  SPURS  OR  2-TONE  3RD  ORDER  IMD? 

TYPE  SPURS  OR  IMD*  SPURS 

DETERMINE  WHAT  KIND  OF  CALCULATION  DESIRED. 

(l)POUT-F(PLO).  (2>P0UT-F(PSIG>.  (3>P0UT«F(PIN> .PIN-PLO-PSIG. 

TYPE  1 t 2 OR  3.»  1 

ENTER  RUN  DATA! 

ENTER  DIODE  911  ALPHA* IS.RS*  40. * 1 . E-12t 10. 

ENTER  INPUT  AND  OUTPUT  IMPEDANCES  ZO.ZL-  50. *50. 

ENTER  FREQUENCIES  FLO*FSIG  (HZ) S 20. .19 
ENTER  POWER  LEVELS  (IN  DBM > * 

PSIGMIN . PSIGMAX * STEPSIG.PLOMIN.PLOMAX  * STEPLO!  -30*-20. 10. 5. 13* 1 
ENTER  NUMBER  OF  SPURS  DESIRED.  (1  TO  10M  A 

ENTER  SPUR  ORDER  M.N.  (M*FLO+N*FSIG> ! l.-l  1.-2  2.-1  2*-2  2*-3  3*-2 

ENTER  SAMPLING  INTERVAL  DELTAT  (SECS)S  l.E-3 

MAX*  NUMBER  OF  ITERATIONS.  MAX.  ALLOWED  ERROR-10. 1 .E-7 

***  END  OF  DATA  INPUT. BEGIN  EXECUTION.  ««* 

«**  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE1 . «** 

»»*  FTWH  TTHF  ftOMATN  OUTPUT  OM  TAPFO  *** 
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* FROM  OOP*  FURNISHED  10  DDC 


EXAMPLE  3 
(b)  Output 


SIMULATION  OF  A BALANCED  MIXER  WITH  A 90.  DEGREE  HYBRID  AND 
IDENTICAL  DIODES 

^CALCULATE  THE  POWER  OUTPUT  OF  SPURIOUS  PRODUCTS 
AS  A FUNCTION  OF  PLO 

DIODEAlt  ALPHA-  40.00  IS-  l.OOE-12  RS-  10.00 
FREQUENCIES:  FLO-  2.00E+01  FSIG-  1 .90E+01  FSIG2-  O.  HZ 

SOURCE  IMPEDANCE-  50.00  LOAD  IMPEDANCE-  50.00 

SPURIOUSM  1.-1H  1.-2H  2.-1H  2.-2X  2.-3X  3.-2) 

RUN  PARAMETERS:  DELTAT-  1.00E-03  NITER-10  ERMAX-  1.0E-07 

*********************************************************************** 


PSIG- 

PLO 

-30.00  DBM 
* ( l.-l)  * 

< 1.-2) 

* 

< 2.-1) 

* 

( 2.-2) 

* < 2.-3) 

* 

< 3.-2) 

* 

5.00 

* 

-41.24 

* 

-91.82 

* 

-48. 5A 

* 

-9A.79 

* -139.51 

* 

-9A.24 

* 

A. 00 

* 

-40. A2 

* 

-97.74 

* 

-49.32 

* 

-93.83 

* -131.53 

* 

-101. 1A 

« 

7.00 

* 

-40.24 

* 

-105.40 

« 

-50. AO 

« 

-92.98 

* -143.00 

* 

-10A.25 

* 

8.00 

* 

-40.00 

* 

-119.18 

« 

-52.34 

* 

-93.04 

* -141.11 

* 

-108. 9A 

» 

9.00 

* 

-39.85 

* 

-111.02 

* 

-54. 5A 

* 

-93.27 

« -128. AO 

* 

-110.11 

* 

10.00 

* 

-39. 7A 

* 

-105. A1 

* 

-57.42 

* 

-94.30 

« -12A.21 

* 

-113.95 

* 

11.00 

* 

-39.71 

* 

-100.87 

* 

-A1.43 

* 

-95.10 

* -12A.14 

* 

-112.83 

* 

12.00 

* 

-39. A9 

* 

-103. A4 

* 

-A7.9A 

* 

-9A.84 

* -14A.27 

* 

-108.98 

« 

13.00 

* 

-39. A9 

* 

-99.33 

* 

-89.58 

* 

-97.47 

* -140.91 

« 

-110.47 

* 

I******************************************************************** 

PSIG- 

PLO 

-20.00  DBM 
* < 1.-1)  * 

( 1.-2) 

* 

( 2.-1) 

* 

< 2.-2) 

* ( 2.-3) 

* 

< 3.-2) 

* 

5.00 

* 

-31.25 

* 

-71.77 

* 

-38.59 

* 

-7A.B8 

* -113.75 

* 

-7A.2A 

* 

A. 00 

* 

-30. A3 

* 

-77.70 

* 

-39.34 

* 

-73.87 

* -108.55 

* 

-81.13 

* 

7.00 

* 

-30.24 

* 

-85.72 

* 

-40. A2 

* 

-73.01 

* -127.74 

* 

-8A.00 

* 

8.00 

* 

-30.00 

* 

-103.01 

* 

-42.35 

* 

-73.07 

* -113.70 

* 

-89.31 

* 

9.00 

* 

-29.85 

* 

-91.55 

* 

-44. 5A 

* 

-73.50 

* -104.31 

* 

-91.02 

* 

10.00 

* 

-29. 7A 

* 

-BA. 04 

* 

-47.43 

« 

-74.34 

* -97. A9 

« 

-93.08 

* 

11.00 

* 

-29.71 

« 

-81. AO 

« 

-51.43 

* 

-75.30 

* -99.12 

* 

-92.07 

* 

12.00 

» 

-29. A9 

* 

-83.75 

* 

-58.04 

* 

-7A.51 

* -121.20 

* 

-91.12 

* 

13.00 

* 

-29. A9 

* 

-79.88 

* 

-81.03 

« 

-77.51 

* -120.39 

* 

-90.51 

* 

*********************************************************************** 


41 


A 


p" 1 f 


I 

I 


1 


AFAL-TR-77-266 


EXAMPLE  4 

CALCULATION  OF  ONE-TONE  IMD  LEVELS  IN  A 180°-HYBRID 
BALANCED  MIXER  AS  A FUNCTION  OF  PSIG 
(a)  Input 


***  NON-LINEAR  CIRCUIT  SIMULATION  *** 

ENTER  OPTIONS  FOR  THIS  RUN} 

ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT. 

(TYPE  1*2*  OR  4 ) ND=2_ 

MIXER  OR  LIMITER?  (TYPE  M OR  LM  M 

90  OR  180  DEGREE  HYBRID?  TYPE  90  OR  1801  180 

IDENTICAL  DIODES?  (YES  OR  NO)}  YES 

TIME  DOMAIN  ONLY?  (YES  OR  N0)1  NO 

CALCULATE  SPURS  OR  2-TONE  3RD  ORDER  IMD? 

TYPE  SPURS  OR  IMD1  SPURS 

DETERMINE  WHAT  KIND  OF  CALCULATION  DESIRED. 

(l)POUT-F(PLO).  ( 2)P0UT=F(PSIG  > • (3)P0UT*»F(PIN)  *PIN=PLO*=PSIG. 

TYPE  1 * 2 OR  3.}_2_ 

ENTER  RUN  DATA1 

ENTER  DIODE  *1}  ALPHA. IS. RS-  40. * 1 .E-12* 10. 

ENTER  INPUT  AND  OUTPUT  IMPEDANCES  ZO.ZL-  50. *50. 

ENTER  FREQUENCIES  FLO.FSIG  (HZ)1  20. *19. 

ENTER  POWER  LEVELS  (IN  DBM) 1 

PLOMIN.PLOMAX* STEPLO* PSIGMIN.PS IGMAX*STEPSIG*  5*10*5*-30*-20»2 
ENTER  NUMBER  OF  SPURS  DESIRED.  (1  TO  10)1  ji. 

ENTER  SPUR  ORDER  M*N.  (M*FLO+N*FSIG) S 1*-1  1.-2  2*-l  2*-2  2*-3  3*-2 

ENTER  SAMPLING  INTERVAL  DELTAT  (SECS)J  l.E-3 

MAX*  NUMBER  OF  ITERATIONS.  MAX.  ALLOWED  ERROR-10* 1 .E-7 

***  END  OF  DATA  INPUT.BEGIN  EXECUTION.  *** 

***  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE1 • *** 

*««  FIND  TIME  DOMAIN  OUTPUT  ON  TAPE2.  *** 


r'-  v. ' mis  PACK  IS  BEST  QUALITY  PRACTICABLE 
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EXAMPLE  4 
(b)  Output 


SIMULATION  OF  A BALANCED  MIXER  WITH  A180.  DEGREE  HYBRID  AND 
IDENTICAL  DIODES 

♦CALCULATE  THE  POWER  OUTPUT  OF  SPURIOUS  PRODUCTS 
AS  A FUNCTION  OF  PSIG 

DIODEfl  ! ALPHA*5  40.00  IS«  1.00E-12  RS«  10.00 
FREQUENCIES  * FLO»  2.00E+01  FSIG=  1.90E+01  FSIG2=  0.  HZ 

SOURCE  IMPEDANCE-  50.00  LOAD  IMPEDANCE-  50.00 

SPURIOUS! ( 1.-1X  1 »-2> ( 2f-l)(  2.-2X  2.-3H  3.-2) 

RUN  PARAMETERS!  DEL TAT-  1.00E-03  NITER-10  ERMAX-  1.0E-07 


****<**********«***************»*«*********************«********«****** 


PLO=  5.00  DBM 


* PSIG 

* 

( l.-l) 

* < 1.-2) 

« 

< 2.-1) 

* ( 2.-2) 

* ( 2.-3) 

» < 3.-2) 

* -30.00 

* 

-43.54 

* -280.09 

* 

-47.89 

* -281.68 

* -128.52 

« -281.57 

« -28.00 

* 

-41.53 

* -277.56 

* 

-45.90 

* -278.85 

* -122.52 

* -279.85 

» -26.00 

* 

-39.53 

* -275.30 

* 

-43.90 

* -276.83 

« -116.40 

« -278.43 

* -24.00 

« 

-37.53 

* -273.82 

* 

-41.90 

* -275.08 

* -110. *0 

» -276.70 

* -22.00 

* 

-35.53 

* -271.88 

* 

-39.91 

* -272.97 

* -104. H6 

* -276.30 

* -20.00 

» 

-33.52 

* -269.81 

* 

-37.92 

* -270.94 

* -98.52 

* -2/5.17 

***«M*«*************«****«****««**«****«*«4t«*»**««****»***»****»«»«*«« 

PLO-  10.00  DBM 

* PSIG 

* 

( l.-l) 

* < 1.-2) 

* 

< 2.-1) 

* < 2.-2) 

» ( 2.-3) 

* < 3.-2) 

* -30.00 

* 

-41.72 

* -279.20 

* 

-50.96 

* -268.07 

* -153.13 

* -273.05 

* -28.00 

« 

-39.72 

* -283.54 

* 

-48.96 

* -266.12 

* -146.52 

* -270.51 

* -26.00 

* 

-37.72 

* -277.68 

* 

-46.96 

* -263.85 

« -137.70 

* -267.93 

« -24.00 

* 

-35.72 

* -271.37 

* 

-44.96 

* -261.94 

* -132.71 

* -266.27 

* >22.00 

* 

-33.72 

* -270.83 

♦ 

-42.96 

« -260.66 

* -126.02 

« -264.97 

* -20.00 

* 

-31.72 

* -272.43 

♦ 

-40.96 

* -258.84 

* -118.71 

« -262.84 

MM****:*************:************************************************** 
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EXAMPLE  5 

CALCULATION  OF  THE  OUTPUT  LEVELS  OF  A 
LIMITER  FED  BY  TWO  SIGNALS 
(a)  Input 


»»t  NON-LINEAR  CIRCUIT  SIMULATION  *** 

ENTER  OPTIONS  FOR  THIS  RUN! 

ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT. 

<TYPE  1 »2»  OR  4)  ND-2 

MIXER  OR  LIMITER?  (TYPE  M OR  L)l  L 

IDENTICAL  DIODES?  (YES  OR  NOM  YES 

TIME  DOMAIN  ONLYT  (YES  OR  NO)>  NO 

CALCULATE  SPURS  OR  2-TONE  3RD  ORDER  IMD? 

TYPE  SPURS  OR  IMDS  SPURS 

DETERMINE  UHAT  KIND  OF  CALCULATION  DESIRED. 

( 1 ) POUT"F ( PLO  > . (2)P0UT-F(PSIG) . (3>P0UT«F(PIN> » PIN«*PLO-PSIG. 
TYPE  1 i2  OR  3.11 

ENTER  RUN  DATA: 

ENTER  DIODE  *i:  ALPHA. IS >RS-  40. . 1 .E-12. 10. 

ENTER  INPUT  AND  OUTPUT  IMPEDANCES  ZO.ZL*  50. .50. 

ENTER  FREQUENCIES  FLO.FSIG  (HZ)J  20. .19, 

ENTER  POWER  LEVELS  (IN  DBM) ( 

PSIGMIN. PSIGMAX.STEPSIG.PLOMIN.PLOMAX.STEPLO:  5. IS. 10. 0.20. 2 

ENTER  NUMBER  OF  SPURS  DESIRED.  (1  TO  10) : 2_ 

ENTER  SPUR  ORDER  M.N.  (M*FLO+N*FSIG) i 1.0. 0.1 

ENTER  SAMPLING  INTERVAL  DELTAT  (SECS)t  l.E-3 

MAX.  NUMBER  OF  ITERATIONS.  MAX.  ALLOWED  ERROR-10. 1 .E-7 

***  END  OF  DATA  INPUT .BEGIN  EXECUTION.  *** 

*««  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE1 • *** 

***  FIND  TIME  DOMAIN  OUTPUT  ON  TAPE2.  *** 
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EXAMPLE  5 
(b)  Output 


SIMULATION  OF  A LIMITER  WITH 
IDENTICAL  DIODES 

♦CALCULATE  THE  POWER  OUTPUT  OF  SPURIOUS  PRODUCTS 
AS  A FUNCTION  OF  PLO 

DIODEtl : ALPHA-  40.00  IS-  1.00E-12  RS-  10.00 
FREQUENCIES i FLO-  2.00E+01  FSIO-  1.90E+01  FSIG2-  0.  HZ 

SOURCE  IMPEDANCE-  50.00  LOAD  IMPEDANCE-  50.00 

SPURIOUS)  ( 1»  OX  Of  1) 

RUN  PARAMETERS:  DELTAT-  1.00E-03  NITER-10  ERMAX-  1.0E-07 


*********************************************************************** 


PSIG- 

PLO 

* 

5.00  DBM 
( 1>  0)  * 

< 

Of  1> 

* 

0.00 

* 

-1.52 

* 

4.07 

* 

2.00 

* 

.33 

* 

3.75 

* 

4.00 

* 

2.13 

* 

3.30 

* 

6.00 

* 

3.87 

* 

2.66 

* 

8.00 

* 

5.55 

* 

1.77 

* 

10.00 

* 

7.11 

* 

.60 

* 

12.00 

* 

8.44 

* 

-.57 

* 

14.00 

* 

9.63 

* 

-1.49 

* 

16.00 

* 

10.79 

* 

-2.25 

* 

18.00 

* 

11.99 

* 

-2.88 

* 

20.00 

* 

13.26 

* 

-3.40 

* 

:«**«****«****««*«****t*******M***t***t**************«**************** 

PSIG- 

PLO 

15.00  DBM 
* < If  0)  * 

< 

Of  1) 

* 

0.00 

* 

-6.92 

* 

10.29 

* 

2.00 

« 

-4.91 

* 

10.27 

* 

4.00 

* 

-2.89 

* 

10.23 

* 

6.00 

* 

-.88 

* 

10.17 

* 

8.00 

* 

1.16 

* 

10.07 

* 

10.00 

* 

3.23 

* 

9.90 

* 

12.00 

* 

5.35 

» 

9.60 

8 

14.00 

« 

7.56 

« 

9.02 

* 

16.00 

* 

9.65 

« 

8.19 

* 

18.00 

* 

11.46 

* 

7.31 

* 

20.00 

* 

12.99 

* 

6.68 

* 

*********************************************************************** 
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EXAMPLE  6 

SIMULATION  OF  A DOUBLE  BALANCED  MIXER 
WITH  INTERMEDIATE  ITERATION  OUTPUT 
(a)  Input 


***  NON-LINEAR  CIRCUIT  SIMULATION  *** 

ENTER  OPTIONS  FOR  THIS  RUN! 

ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT. 

(TYPE  1*2*  OR  4)  ND-4 

IDENTICAL  DIODES?  (YES  OR  NO) l YES 

TIME  DOMAIN  ONLY?  (YES  OR  NO)!  YES 

PRINTOUT  OPTIONS? 

INTERMEDIATE  ITERATION  OUTPUT?  (YES  OR  NO)!  YES 
ENTER  RUN  DATA! 

ENTER  DIODE  fit  ALPHA* IS*RS=  40. * 1 . E-12. 10 . 

ENTER  INPUT  AND  OUTPUT  IMPEDANCES  ZO.ZL=  50.* 50. 

ENTER  FREQUENCIES  FLO*FSIG  (HZ):  20. *19. 

ENTER  POWER  LEVELS  (IN  DBM)! 

PLO.PSIG-  10.-10 

ENTER  SAMPLING  WINDOW  TMAX  (SECS)-  l.E-2 

ENTER  SAMPLING  INTERVAL  DELTAT  (SECS)!  2.E-3 

MAX*  NUMBER  OF  ITERATIONS*  MAX.  ALLOWED  ERROR-10, 1 ,E-7 

***  END  OF  DATA  INPUT  * BEG IN  EXECUTION.  *** 

***  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE1 . *** 
tit  FIND  TIME  DOMAIN  OUTPUT  ON  TAPE2.  *** 
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EXAMPLE  6 
(b)  Output 


SIMULATION  OF  A DOUBLE  BALANCED  MIXER  WITH 
IDENTICAL  DIODES 

DIODEfl * ALPHA=  40.00  IS"  1.00E-12  RS=  10.00 
FREQUENCIES!  FLO"  2.00E+01  FSIG=  1.90E+01  FSIG2=  0.  HZ 

SOURCE  IMPEDANCE"  50.00  LOAD  IMPEDANCE"  50.00 

RUN  PARAMETERS!  DELTAT"  2.00E-03  NITER"10  ERMAX=  1.0E-07 


**M*:M*********«*************************«**************************** 


PSIG=  -10. 
T VD( 1 > 

00  DBM  PLO=  10.00  DBM 

VB<2>  VD<3)  VD<4) 

ERROR 

ITER 

IDT 

0. 

0. 

0.  0.  0. 

0. 

1 

1 

0. 

0. 

0.  0.  0. 

0. 

2 

1 

2.000E-03 

2.250E-01 

2 . 723E-01 -2 . 250E-01 -2 . 723E-01 

2.49AE-01 

1 

1 

2.000E-03 

2.250E-01 

2. 723E-01-2. 250E-01-2. 723E-01 

5. 5A0E-1 1 

2 

1 

4.000E-03 

4 . 358E-01 

5 . 277E-01 -4 . 358E-01-5 . 277E-01 

2. 192E-01 

1 

1 

4 • 000E-03 

4 . 470E-01 

5 • 0AAE-01-4 .213E-01-5.351E-01 

8 . 338E-04 

2 

1 

4.000E-03 

4.53AE-01 

4.952E-01-4. 13AE-01-5.395E-01 

2.495E-04 

3 

1 

4.000E-03 

4 . 547E-01 

4 . 932E-01 -4 . 122E-01-5 . 402E-01 

8. 003E-0A 

4 

1 

4.000E-03 

4.547E-01 

4.931E-01-4.121E-01-5. 402E-01 

4 . 903E-09 

5 

1 

A.000E-03 
A.000E-03 
A.000E-03 
A.000E-03 
A.000E-03 
A . 000E-03 
6.000E-03 
A.000E-03 
A.000E-03 


A. 7A5E-01 
A. 51AE-01 
A.2A8E-01 
A . 027E-01 
5.802E-01 
5. A20E-01 
5.517E-01 
5.492E-01 
5. 491E-01 


5.981E-01-5. 143E-01-7.873E-01 
5. 758E-01-4 • 983E-01-7 . A9AE-01 
5 . 58AE-01-4 . 859E-01 -7 . 520E-01 
5.512E-01-4.803E-01-7.351E-01 
5.51AE-01-4.801E-01-7. 195E-01 
5 . 528E-01-4 . 80AE-01-7 . 0A8E-01 
5 . 535E-01-4 . 809E-01-A. 997E-01 
5 . 536E-01-4 . B09E-01-6 . 980E-01 
5 . 53AE-01-4 . 809E-01 -6 . 979E-01 


1.317E-01 
1.695E-03 
1.3A7E-03 
9.567E-04 
7 • 500E-04 
4 • 931E-04 
1 . 552E-04 
9.220E-06 
2.096E-08 


1 

2 

3 

4 

5 
A 

7 

8 
9 


1 

1 

1 

1 

1 

1 

1 

1 

1 


8.000E-03  5 . 7A2E-01 
8.000E-03  5. A93E-01 
8.000E-03  5. A82E-01 
8 . 000E-03  5. A82E-01 


5. 774E-01-5. 1A8E-01-7.89AE-01 
5.71AE-01-5. 127E-01-7.847E-01 
5.709E-01-5. 122E-01-7.839E-01 
5.709E-01-5. 122E-01-7.839E-01 


1 . 099E-02 
1 . 215E-04 
2.435E-0A 
B.219E-10 


1 1 
2 1 

3 1 

4 1 


I*******************************************************:*************** 
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6.  PROGRAM  LISTING 

The  comments  on  the  enclosed  listing  make  the  program  MIXER  self 
explanatory.  The  equivalence  between  various  symbols  in  the  program 
and  in  the  test  follows: 


PROGRAM 

AJ 

AJINV 


TEXT 

J 

J" 


FV 


F 


VD 


v 
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PROGRAM  M I XER ( I NPUT  , OUTPUT  » TAPE 1 , TAPE2 ) 

C ****************************************************************** 
C * THIS  PROGRAM  CALCULATES  THE  PERFORMANCE  OF  NON-LINEAR  CIRCUITS 
C « CONTAINING  SEMICONDUCTOR  DI0DE8.  (SINGLE  ENDED , BALANCED » 

C * AND  DOUBLE  BALANCED  MIXERS*  AND  BACK  TO  BACK  DIODE  LIMITER8). 

C * THE  VARIOUS  OPTIONS  AREt 
C * IDENTICAL  OR  NON-IDENTICAL  DIODES. 

C * SOLUTION  IN  THE  TIME  DOMAIN  ONLY. 

C * CALCULATION  OF  ONE  TONE  IMD '8  OR  TWO  TONE  THIRD  ORDER  IMD. 

C * ONE  OF  THE  FOLLOWING  POWER  COMBINATIONS! 

C * A) CONSTANT  SIGNAL  POWER I POUT  AS  A FUNCTION  OF  LO  POWER. 

C * B>  CONSTANT  LO  POWER!  POUT  AS  A FUNCTION  OF  SIGNAL  POWER. 

C * C>  SIGNAL  POWER-LO  POWER. 

C ****************************************************************** 

c 

IMPLICIT  LOGICAL(L) 

REAL  IS(4>,ALPHA(4>,RS(4),E(4),E1(4>,EP(4>,VD(4>,ID(4> 

DIMENSION  AJ(4,4),AJINV(4,4),WK(4) t HPROD< 10) ,NPROD( 10) ,OMEG< 10)  , 
1P0UT  (10). VS INC 10) , VCOS( 10  > , CVD( 4,4)  , CE(4,4),CVIN(4,4), 

1VIN(4) f UFINAL(4> .FV<4> 

COMMON/Z/ZO 

PRINT  1 a 

LSAME-LTDOM-LTHIRD-LITER~LIMIT~. FALSE.  \ 

C I 

C ******************************************************************* 
C * BEGIN  OPTIONS  INPUT  * 

C ******************************************************************* 

PRINT*.’  ENTER  OPTIONS  FOR  THIS  RUN!’ 

PRINT*.’  ENTER  NUMBER  OF  DIODES  FOR  SIMULATED  CIRCUIT.* 

PRINT*.’  (TYPE  1.2.  OR  4)  ND-* 

READt.ND 

IF(ND.EQ.4)G0  TO  90 
LSAME-ND.EG. 1 
IF (ND.EQ.  1 )G0  TO  95 

PRINT*.’  MIXER  OR  LIMITER?  ( TYPE  "M  OR  L>!  • 

READ  2.ANS 
LIMIT-ANS.EQ.1HL 

IF< .NOT .LIMIT) PRINT*. * 90  OR  180  DEGREE  HYBRID?  TYPE  90  OR  1801’ 
IF ( .NOT .LIMIT)READ*.DEG 
90  PRINT*,’  IDENTICAL  DIODES?  (YES  OR  N0>1  • 

READ  2.ANS 
LSAME~ANS.EQ.1HY 

95  PRINT*,’  TIME  DOMAIN  ONLY?  (YES  OR  NO) I • 

READ  2.ANS 
LTD0M~ANS.E0.1HY 
IF(LTD0M)00  TO  100 

PRINT*,*  CALCULATE  SPURS  OR  2-TONE  3RD  ORDER  IMD?  * 

PRINT*,’  TYPE  SPURS  OR  IMD1  • 

READ  2.ANS 
LTHIRD-AN8.EQ. 1HI 

PRINT*,’  DETERMINE  WHAT  KIND  OF  CALCULATION  DESIRED.* 


quality  mcncAMi 
ED  TO  DOC  - 
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PRINT*.*  < 1 ) POUT -F ( PLO  > • <2>P0UT-F!P8IG> . • . 

1*  <3)P0UT-F(PIN) .PIN-PL0-P8IG. * 

PRINT***  TYPE  1 .2  OR  3.1* 

READ** I FLAG 
GO  TO  110 

100  PRINT***  PRINTOUT  OPTIONS!  * 

PRINT***  INTERNEDIATE  ITERATION  OUTPUTT  (YES  OR  NO)!  * 

READ  2*ANS 

I_ITER-ANS.EQ.1HY 

IFLAG-1 


« BEGIN  DATA  INPUT  * 

110  PRINT** * ENTER  RUN  DATA!* 

ND1-ND 

IF (LSAME)NDl-l 
DO  115  1-1 *ND1 

PRINT***  ENTER  DIODE  •*.I**!  ALPHA. IS*RS-  * 

115  READ*  * ALPHA  <I)*IS(I)*RS(I) 

IF< .NOT «LSAHE>GO  TO  140 
DO  120  1-2*4 
ALPHA! I ) -ALPHA! 1 ) 

IS!I)-IS!1) 

120  RS! I )-R8! 1 ) 

140  PRINT*.*  ENTER  INPUT  AND  OUTPUT  IMPEDANCES  70. ZL-  * 

READ«*ZO*ZL 

PRINT*. * ENTER  FREQUENCIES  FLO.FSIG  !HZ)!  * 

READ* * FLO* FSIG 

IF !LTHIRD> PRINT* * * ENTER  2ND  TONE  FREQUENCY  FSI02!  * 
IF!LTHIRD>READ*.FSI02 

PRINT*.*  ENTER  POWER  LEVELS  !IN  DBM)!* 

IF !LTDOM)GO  TO  152 
IF! IFLAG.EQ. 3)00  TO  144 
IF! IFLAQ.EQ. 1 ) PRINT*. 

1 ' PSIGMIN.PSIGMAX . STEPS IG * PLOM IN. PLOMAX . STEPLO! * 
IF!IFLAG.EQ.2)PRINT** 

1*  PLOMIN. PLOMAX. STEPLO* PSIGMIN. PSIGMAX.STEPSIGt * 

READ* . PCMIN . PCMAX . PCSTEP . PMIN, PMAX . PSTEP 
GO  TO  148 

144  PRINT***  PMIN.PMAX. PSTEP-  * 

READ*. PMIN. PMAX. PSTEP 

148  IF!LTHIRD)00  TO  156 

PRINT***  ENTER  NUMBER  OF  SPURS  DESIRED.  (1  TO  10)1  * 

READ»*NPRODS 

PRINT*.*  ENTER  SPUR  ORDER  M.N.  (M*FL0+N*F8I0) 1 * 

READ*  * ! MPROD ! I ) * NPROD ! I > * I - 1 * NPRODS  > 

GO  TO  156 

152  PRINT*.*  PLO.PSIG-  * 
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READ* » PMI N » PCM I N 
PCMAX-PCMIN 
PMAX-PMIN 
PCSTEP-PSTEP-O. 

PRINT*# * ENTER  SAMPLING  WINDOW  TMAX  ( SEC8 ) ■ • 

READ*»TMAX 

136  PRINT*# * ENTER  SAMPLING  INTERVAL  DELTAT  (SECS>1  • 

READ*# DEL TAT 

PRINT*.*  MAX#  NUMBER  OF  ITERATIONS#  MAX.  ALLOWED  ERROR-  * 

READ*  # NITER  # ERMAX 

PRINT*. * ***  END  OF  DATA  INPUT.BEOIN  EXECUTION.  ***• 

PRINT*.*  ***  FIND  FREQUENCY  DOMAIN  OUTPUT  ON  TAPE1 • **** 

PRINT*. * ***  FIND  TIME  DOMAIN  OUTPUT  ON  TAPE2.  **** 

C 

C ******************************************************************* 
C * BEGIN  TITLE  PRINTOUT  » 

C ******************************************************************* 

I TAPE-1 

IF<LTDOM> ITAPE-2 

IF! ND.EQ.l) WRITE ( ITAPE.O) 

IF! ND.EQ.4) WRITE! I TAPE# 10) 

IF (LIMIT ) WRITE ( I TAPE  #12) 

IF(ND.EQ.2. AND. .NOT • LIMIT )WRITE( I TAPE  > 14>DE0 
IF (L8AME.AND.ND.NE. 1 >WRITE( I TAPE# 14) 

IF < .NOT .LSAME) WRITE  < I TAPE  #18) 

IF(LTD0M)00  TO  148 
IF< IFLAO.EQ. 3 >G0  TO  166 
PV-4H  PLO 

IF( IFLAO.EQ. 2)PV-4HPSI0 
IF(LTHIRD)GO  TO  162 
WRITE ( 1 #24)PV 
00  TO  168 
162  WRITE< 1 #26>PV 
00  TO  168 
146  WRITE ( 1 #28) 

168  CONTINUE 

DO  170  I-l.NDl 

170  WRITE< I TAPE #30) I # ALPHA ( I >#IS(I)»RS(I) 

WRITE ( I TAPE  »32)FL0»FSI0»FSIG2 
WRITE! ITAPE #33)Z0#ZL 

IF< .NOT • < L THIRD. OR. LTDOM) ) WRITE! 1 #34)NPR0DS» 

1 (MPROD! I > » NPROD ! 1 > # 1-1 .NPRODS  > 

WR I TE ( I TAPE # 36 > DEL T A T # N I TER # ERMAX 

******************************************************************* 
* BEGIN  INITIALIZATION  * 

IF (LTDOM) NPRODS- 1 
IF(LTMIRD)NPR0DS-2 
PI-AC0S(-1 . > 
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RT2-SQRT(2. ) 

0MSIG-2*PI*FSIG 

0ML0-2.*PI*FL0 

NP-NPC-1 

IF  < PSTEP . NE . 0 . > NP-1  + < PMAX-PHIN ) /PSTEP 
IF ( PCSTEP . NE . 0 . )NPC-1+  < PCMAX-PCMIN ) /PCSTEP 
IF(LTDOH)GO  TO  230 
TMAX-1 . /ABS  < FLO-FSIG ) 

IF<LTHIRD)GO  TO  200 
DO  180  I-l f NPRODS 

1 80  OMEG ( I ) -ABS ( MPROD ( I ) *OML O+NPROD < I > *OMSIG ) 

GO  TO  220 

200  TNAX-AHAX1 ( TMAX r 1 • /ABS(FSIG-FSIG2> > 

0MSIG2-2. *PI*FSIG2 

OHEG ( 1 > -ABS  < 2 . *OMS IG-ONS IG2-0ML0 > 

OMEG ( 2 ) -ABS  < OMLO-OMSIG ) 

220  IF< IFLAG . NE .3)00  TO  230 
APV-4H  PIN 
URITE( 1 >38) 

IF ( . NOT . LTH I RD ) UR I TE  < 1 f 40 ) APV  t NPRODS , ( MPROD (I)r NPROD (I)r 
11-lr NPRODS) 

IF (LTHIRD)URITE( 1 r42)APV 
C 

C ******************************************************************* 
C * INITIALIZE  CIRCUIT  MATRICES  * 

C ******************************************************************* 

230  DO  240  I-lr4 
DO  235  J-lr4 

235  CVD( I r J >— CE ( I r J )— CVIN ( I r J ) — 0 ♦ 

240  CONTINUE 

IF(ND.EQ.4)G0  TO  250 
IF(ND.E0.2)G0  TO  242 

**«*****»*»*«»*«*»«*«* 

* SINGLE  ENDED  MIXER  * 

********************** 

CVD< 1 r 1 )-l . 

CE ( 1 r 1 > -R8 ( 1 > +Z0 
CVIN(lrl)  — 1. 

GO  TO  260 

****************** 

* BALANCED  MIXER  « 

«»***«********»*«« 

242  CVD( 1 r 1 )>CVD(2r 2)-l • 

IF (LIMIT )G0  TO  245 
CE ( 1 r 1 ) -ZL +ZO+RS ( 1 > 

CE ( 2 r 2 > -ZL+ ZO+RS ( 2 > 

CE ( 1 r 2 ) -CE ( 2 r 1 ) — ZL 
CVIN(lrl)— 1. 
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CVIN(2f2>-1. 

00  TO  260 

**»****»*«« 

* LIMITER  * 

««**»*****» 

245  DIV-ZL/<ZL+Z0) 
Z1-DIV4Z0 
CE( 1 r 1 >"RS< 1 )+Zl 
CE  < 2 * 2 ) “RS < 2 ) +Z1 
CE(1»2>-CE<2»1>  — Z1 
CVINU'l)  — DIV 
CVIN(2>2)-DIV 
00  TO  260 


***»****»t**»«**********« 

* DOUBLE  BALANCED  HIXER  * 

»«*»**«*»»«*******«**«*** 

250  CVD< 1 » 1 )«CVD< 1 r 2)«CVD( 2» 2 )“CVD< 2» 3 ) “CVD< 3 » 2> "CVD( 4 » 1 )■ 
lCVD<4r2)-CVD<4r3)-CVD<4.4>-l. 

CE ( 1 > 1 ) «RS < 1 ) +Z0/2 . 

CE ( 1 r 2 > *CE  < 2 » 2 ) *RS  < 2 > +Z0/2 • 
CE<1f3)-CE<1f4)-CE<2»1>-CE<2.4>— ZO/2. 

CE  < 2 > 3 ) «RS < 3 > +Z0/2 . 

CE(3»1>»CE(3f3)— (Z0/2.+ZL) 

CE  < 3 > 2 ) »RS < 2 ) +Z0+ZL 
CE(3t4)-ZL 
DO  255  I-lt4 
255  CE<4»I)-R8(I) 

CVINU'l)  — 1. 

CVIN<2'2>  — 1. 

CVIN(3'1)-CVIN<3'2>— .5 
260  CONTINUE 


PARAMETER  LOOP 


DO  1000  KCBTEP-l'NPC 

PC0N8T  *PCMIN+  < KC8TEP- 1 )*PC8TEP 

VC0N8T  «F  < PCONST  > 


* VARIABLE  POWER  LOOP  * 


DO  900  KSTEP-l.NP 
T-0. 

DO  280  I-1'NPRODS 
280  VC08<I)-VSIN<I>-0. 
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DO  290  I-1fND 

290  VIN(I)-VFINAL(I)-VD(I>-0. 

PVAR-PMIN+(K8TEP-1 >*PSTEP 
WAR-F(PVAR) 

IF(IFLAG.EQ.3)G0  TO  320 
IF(IFLAG.EQ.2)G0  TO  300 
VSIQ-VPONST 
VLO-WAR 
P8I0-PC0NST 
PLO-PVAR 
APC-3HPSIG- 
APV-4H  PLO 
00  TO  340 
300  VS  10- WAR 
VLO-VCONST 
PLO-PCONST 
PSIO-PVAR 
APC-5H  PLO- 
APV-4HPSI0 
00  TO  340 
320  VS  10- WAR 
VLO-VSIG 
PSIG-PLO-PVAR 
340  VSI02-0. 

IF  <LTHIRD)VSIG2-VSIG 

IF (LTDOM. OR. <KSTEP .GT . 1 ) .OR. ( IFLAG .EQ. 3 ) >00  TO  345 

UR I TE  < 1 * 44 ) APC * PCONST 

IF  < .NOT .LTHIRD)WRITE(1f40> 

I APV r WPRODS r ( HPROD (Dr NPROD ( I ) r I -1 r NPRODS ) 
IF<LTHIRD)URITE( 1 »42)APV 
345  IF(LTD0M)URITE(2f46)PSI0fPL0 

IF (LTDOM . AND . .NOT .LITER)URITE< 2f 48) 

ITER-IDT-1 

IF (LITER) WRITE < 2f49) 

****************** 

* TIME  STEP  LOOP  * 

****************** 

350  CONTINUE 

CALCULATE  INPUT  VOLTAGES  AT  T 

IF  (ND.EQ  .4  >00  TO  370 
IF (LIMIT. OR. ND.EQ. 1)00  TO  360 

BALANCED  MIXER 

IF ( DEO. EQ. 180. )G0  TO  354 

VIN( 1 )-RT2/2.*( VLO*SIN(OMLO*T) +VSIG*COS(OMSIG*T>+ 

1 VS IG2*C0S ( 0 MS 102*  T > ) 
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WXS  BJ$T  QUALITY  PRACTICABLE 

,|8QK  tXiPxTUR|lllSHED  TO  DDC  


VIN( 2)“RT2/2.*< VLO*COS<  OMLO*T ) +USIG*SIN(OHSIG*T )+ 

1 VSIG2*SIN(0MSIG2*T> ) 

IF< T.GT.O. >GOTO  380 
C 

C INITIAL  GUESS  AT  T-0 
C 

VD<1)«VIN<1> 

VD<2>— VIN<2> 

VMAX-1/ALPHA< 1 > <ALOG<ABS< VIN< 1 ) )/(CE< 1 t 1 >*IS( 1 > ) > 

IF  < VD  < 1 ) • LT • UMAX > GO  T0380 
VD<  1 >*VMAX 
GO  TO  380 

354  VIN1-RT2/2.*VL0*SIN<0ML0*T> 

VIN2-RT2/2.*<VSIG*SIN<0MSIG*T>+VSIG2*SIN<0MSIG2*T> ) 

VIN<1>— VIN1+VIN2 
VIN(2)-VIN1+VIN2 
GO  TO  380 
C 

C LIMITER  AND  SINGLE  ENDED  MIXER 

C 

360  VIN< 1 >"VIN<  2 )“VLO*SIN<OMLO*T ) +VSIG*SIN<  OMSIG*T ) + 
1V8IG2<SIN<0MSIG2<T > 

GO  TO  380 
C 

C DOUBLE  BALANCED  MIXER 

C 

370  VIN< 1 )"VLO*SIN<OMLO*T ) 

VIN<2)«VSIG*SIN<0MSIG*T)+VSIG2*SIN<0MSIG2*T> 

380  CONTINUE 
C 

C ******************************************************************** 
C 8 ITERATION  LOOP  * 

C ******************************************************************** 

c 

DO  390  I-ltND 
El<I>-0. 

IF< ALPHA < I )*VD< I ) .GT . 100. )G0  TO  420 

IF< ALPHA< I >*VD< I ) .GT. -100. )E1<I)«EXP( ALPHA  < I > * VD  < I ) ) 
E(I)-I3<I)*(E1<I)-1.) 

390  EP< I >«IS< I >«ALPHA( I )SE1 ( I ) 

C 

C CALCULATE  FV  AND  JACOBIAN 

C 

DO  400  I-1'ND 
FV< I >“0. 

DO  393  J-1'ND 

FV<  I )«FV<  I )+CVD(  I » J)*VD<  J) +CE< I f J)*E< J) +CVIN<  IrJ)*VIN<J) 

395  AJ< I » J>“CVD< I » J>'fCE< I » J ) *EP( J) 

400  CONTINUE 
C 
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THIS  PAGE- 1$  BEST  QUALITY  PRACTICABLE 
JHOH  COPY 'FUtotl SHED  TO  DDC  


c 

c 

c 


CALCULATE  INVERSE  OF  JACOBIAN 

CALL  LINV1F< AJ>NDt4r AJINVrO»WK>IER) 

THE  MATRIX  AJ  IS  SINGULAR  WHEN  IER-129 

IFdER.EQ.  129)00  TO  440 

ERROR-O. 


C 

C CALCULATE  DIODE  VOLTAGES 

C 

DO  410  I-1>ND 
VDN-O. 

DO  40S  J-l.ND 

405  VDN-VDN+AJINV(I.J)*FV< J) 

ERR0R-ERR0R+VDNM2 
VDN-VD(I)-VDN 
410  VD( I >-VDN 

IF  <LITER)WRITE<2'60)T  .VD. ERROR. ITER. IDT 
IF(ERROR.LT.ERMAX. AND. ITER. GT.1JG0  TO  440 
ITER-ITER+1 

IF( ITER.GT .NITER) GO  TO  420 
GO  TO  360 
C 

C IF  SOLUTION  DOES  NOT  CONVERGE  REDUCE  TINE  STEP 

C 

420  ITER-1 

IF < T . EQ . 0 • >G0  TO  440 
T-TLAST+DELTAT/<  2**IDT ) 

IF<T.LE.T1  >T-Tl+DELTAT/<2«< IDT+1 ) > 

DO  430  1-1 »ND 
430  VD< I )-VFINAL< I > 

IDT-IDT+1 

IF< IDT .GT .4)00  TO  440 
GO  TO  350 

440  WRITE( ITAPE.45)T 
PRINT  45.  T 
STOP 

440  IF (LITER >WRITE< 2.50) 

C 

C CALCULATE  DIODE  CURRENTS  AND  LOAD  VOLTAGE 

C 

DO  4B0  I-l.ND 
ID<I>— 18(1) 

480  IF(ALPHA(I)*VD<I).GT.-100.) 

1 1 D < I ) - IS  < I ) * ( EXP  < ALPHA  < I > *VD ( I ) ) - 1 . > 

VLOAD-C ID< 1 )-ID(2> )*ZL 

IF(ND.EQ. 1 )VLOAD-ID< 1 )*RS< 1 >+VD< 1 )-VIN< 1 )/2. 

IF (LIMIT  >VLOAD-VD< 1 ) + ID< 1 >*RS( 1 ) 

IF (ND.EO. 4 ) VLOAD-ZL*< ID(4)4ID(2>-ID(1)-ID(3)> 

IF (LTDOM.AND. . NOT. LITER)MRITE<2»40)T»VD*VL0AD» ITER. IDT 
IF < IDT.GT . 1 >G0  TO  520 


j 
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YWfcRMEriS  BEST  QUALITY  PRAOTIOASM 
tXM  OOPY 'FURNISHED  TO  DDC  


C 

c 

c 


TLAST-T 

IF<LTDOM)GO  TO  32 0 

CALCULATE  FOURIER  COMPONENTS  AND  POWER  OUTPUT 
DO  500  MN-lr NPRODS 

VCOS ( MN ) -VCOS ( MN  > +VLOADBCOS ( OMEO < MN  > BT ) 

500  VSIN<MN)-VSIN<MN)+VLOADSSIN<OMEG<MN>BT> 

320  Tl-T 

DO  530  I-lrND 
530  UFINAL ( I )-VD( I ) 

T-TLA8T+DELTAT 

ITER-XDT-1 

IF<T.LT.THAX-DELTAT/10.)G0  TO  330 
ADT- < 2 . BDEL  TAT/TMAX > **2 
DO  560  HN-lrNPRODS 
IF(LTD0M>00  TO  900 

P0UT(HN)-10.*AL0010(ADT*<UC08(MN)«*2'FUSIN(HN>**2)/(2.*ZL)*1.E3) 
IF<0ME0<MN) . EQ.O. )P0UT(MN)-P0UT(MN)+3. 

560  CONTINUE 

WRI TE  < 1 r 70  ) PUAR  r NPRODS  r < POUT  ( MN  > r MN- 1 r NPRODS  > 

900  CONTINUE 
1000  CONTINUE 

WRITE < ITAPEr 72) 

IF  < .NOT  <LTDOM)WRITE< 1 r 76 ) 

STOP 


C 

C 


* FORMAT  STATEMENTS 

***»»**«**t«*M****M**««tt«*««t***«*************»****«*«**»****1 

1 FORMAT < * ***  NON-LINEAR  CIRCUIT  SIMULATION  «***r/> 

2 FORMAT <A1) 

8 FORMAT <1H1.3X. 'SIMULATION  OF  A SINGLE  ENDED  MIXER*) 

10  FORMAT < 1 HI r5Xr ‘SIMULATION  OF  A DOUBLE  BALANCED  MIXER  WITH*) 

12  FORHATdHlrSXr ‘SIMULATION  OF  A LIMITER  WITH*) 

14  FORHATdHlrSXr ‘SIMULATION  OF  A BALANCED  MIXER  WITH  A*rF4.0r 
1*  DEGREE  HYBRID  AND*) 

16  FORMAT <6Xr ‘IDENTICAL  DIODES* r/) 

18  F0RMAT(6Xr ‘NON-IDENTICAL  DIODES*  r/> 

24  FORMAT (SXr'BCALCUL ATE  THE  POWER  OUTPUT  OF  SPURIOUS  PRODUCTS* r 
l/5Xr  * AS  A FUNCTION  OF  *rA4r/> 

26  FORMAT < SXr'BCALCUL ATE  THE  POWER  OUTPUT  OF  2-TONE  3RD  ORDER* r 
1*  IHD  PRODUCTS* r/6Xr 'AS  A FUNCTION  OF  *rA4r/> 

28  FORMAT (5Xr*BCALCUL ATE  THE  POWER  OUTPUT  OF  SPURIOUS  PRODUCTS *r 
1/SXr*  AS  A FUNCTION  OF  PIN.PIN-P8IG-PL0. * r/) 

30  FORMAT (3Xr ‘DIODE#*  r II r ' t ALPHA-* rF6.2r3Xr ' IS-* r lPE9.2r 
13Xr*RS-*rOPF7.2) 

32  FORMAT (5Xr ‘FREQUENCIES!  FLO-* r lPE9.2r 2Xr ‘F8IG-* r lPE9.2r 
12Xt ‘F8IG2-*  r 1PE9. 2r • HZ*) 

33  FORMAT <5Xr* SOURCE  IMPEDANCE- • ,F6. 2 r lOXr ‘LOAD  IMPEDANCE-* r 
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25 18  BBST  QUALIrY  mcncABLa 

OQPy  jfURHISHgD  no  jjjjq  


1F6.2) 

34  F0RMAT(5X, 'SPURIOUS! • »•<  * < ' , 12, * . ' « 12, * > * > > 

36  FORMAT (SXt* RUN  PARAMETERS t DELTAT-' , 1PE9.2,3X, 'NITER-' , 
112, SX,'  ERMAX-',lPE8.1,/> 

38  FORMAT </,lX,120<*8* )//5X»  *PIN*PL0"P8IG' > 

40  FORMAT < • ** ,3X, A4,2X, '*• ,-< IX, * < * , 12, * , * , 12, * > *')) 

42  FORMAT ( • *• ,3X,A4,2X, **  IHD  * IF  **) 

44  FORMAT </, IX, 120 < '*' >//3X,A5,F7*2, * DBM') 

46  FORMAT (/, IX, 120 ( •*'>//10X»'P8I0"'»F7.2*'  DBM', 
13X»'PL0-',F7.2,'  DBM') 

48  FORMAT <3X, 'T',T12,'VDC1> • ,T22, *WD<2) ' , T32, 'VD<3>  * » 

1T42, *VD<4> ' ,T32, 'VLOAD* ,T62, 'ITER' , IX, ■ IDT' ) 

49  FORMAT < 3X, *T* ,T12, *WD< 1 ) * , T22, *VD<2) ' ,T32, *WD<3) ' , 

1T42, 'VD<4> ' , T32, 'ERROR* ,T62, 'ITER' , IX, 'IDT*  > 

50  FORMAT</> 

60  FORMAT (6( 1PE10, 3) ,2X,2I3) 

65  FORMAT < ' THE  SOLUTION  DOES  NOT  CONVERGE  AT  T-'»1PE15.4> 
70  FORMAT  < * *• ,F7.2,2X, ,-< 1X,F7.2, ' *•>> 

72  FORMAT (/,1X,120('8*>> 

74  FORMAT  < 1H1 ) 

END 


FUNCTION  F<P > 

COMMON/Z/ZO 

F-10.**(P/20)*SQRT(8.*Z0/1000. > 

RETURN 

END 
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